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1 Introduction 



1.1 Physics Background and 6is 

Since neutrinos have peculiar properties such as very smaU masses, very much 
different mixing patterns from quark sector and possibility of being Majorana 
particle, study of neutrinos is very important to understand the nature of the 
elementary particles and to construct its unified theory in the future. Because 
the neutrinos have finite masses and mixing, they perform neutrino oscillation 
if their masses are not identical. For two flavor case, the probability that a 
particular flavor neutrino with energy E,^ to remain as after traveling a 
distance L becomes, 

P.,-.., = l-sin2 20sin"$, (1) 
Am^L 



due to the neutrino oscillation, where Am^ is the difference of the squared 
masses between the two mass eigen states; — ml, and 9 is the mixing angle 
between mass eigenstate and flavor eigenstate. Through neutrino oscillation, it 
is possible to access to a very small mass range where direct measurements are 
difficult to reach and to measure the mixing angle. 

For 3 flavor case, the flavor eigenstates; Ve, and Vr, and mass eigenstates; 
vi, V2 and z^s, are mixed by 3 x 3 Maki-Nakagawa-Sakata (MNS) matrix jlj. 

Uel Ue2 Ue3 

(7^1 U,,3 I I ^^2 I = Umns I 1^2 I (2) 

UtI Ur2 Ur3 

In this case, three flavor oscillation formula becomes, 
= Sc.f3-4.J2 ^iU*,Ufs,U^,U*p^) sin2 ± 2 ^ 3([/*,[/^,;7„,[/^*^.) sin 2$^ 

i>j i>3 

(3) 

where a, (3 are indices of flavors and i,j are indices of masses with being the 
mass of Vi. <i>y = ^^'^ represent oscillation phase caused by the difference of 
mass square, Am^^- = m'j — mf. 







V, 
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The MNS matrix can be expressed by 4 parameters; 



/I \ / ci3 siae'^' \ / \ 

UmNS = C23 S23 1 -Si2 C12 

V -S23 C23 / V -«i3e"*''' ci3 y V 1/ 



where and Sy express cos 0y and sin 6ij , respectively. 

Experimental studies of the neutrino oscillations have been progressing very 
rapidly in these days. The first firm evidence of the neutrino oscillation was 
discovered in the disappearance of muon type atmospheric neutrinos by Su- 
perKamiokande (SK) group in 1998 0. SK group has been taking more data 
and recently observed spectrum distortion [3], which is consistent with oscilla- 
tion pattern. K2K group confirmed the result 0] by observing the disappearance 
of t'p produced by the KEK-PS accelerator, after traveling 250 km. These phe- 
nomena are naturally explained by the i'^ — > Vr oscillation since there is no 
significant enhancement in i/^ spectrum. The measured oscillation parameters 
are Am23 ~ 2.4 x 10~^eV^ and sin^ 2023 ~ 1. As for the electron type neutrinos 
there have been indications of the neutrino oscillation in the deficit of the solar 
neutrinos for a long time [S] . A transformation of the solar v^. to other type neu- 
trinos was identified by combining SNO data and SK data in 2002 PlITl. Finally, 
the KamLAND group observed a large deficit in reactor De in 2002 |8j and then 
spectral distortion in 2005 Combined with solar neutrino data, the oscilla- 
tion parameters are measured as Am^j ~ '''•9 ^ 10~^eT^^ and tan^ 9i2 « 0.40. 
With these measurements, measurable 2 mass combinations and 2 mixing an- 
gles out of 3 are already measured. The sign of Am^j is determined to be > 
by MSW matter effect of solar neutrinos. However, the sign of Am23 is not 
determined yet. Finite ^13 and Si have not been measured yet. The most strin- 
gent upper limit of ^13 was measured by CHOOZ group using reactor D^s to be 
sin^ 26*13 < 0.15 at Am^ = 2.5 x 10~^eV^ [EI. As for Si parameter, there is 
no information by now. By combining measured parameters, the magnitudes of 
MNS matrix elements are determined to be roughly. 



and the relation of masses and flavor components of mass eigenstates are shown 
in Fig. n The elements of the MNS matrix are very much different from those 
of the Cabbibo-Kobayashi-Maskawa Matrix, which is nearly diagonal. The Ues 
element (and thus sin^is) is small compared to other elements. This peculiar 
feature has to be explained by the unified theory of elementary particles and the 
smallness of the ^13 may play a key role when constructing the unified theory. 




(4) 
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Normal Hieraicliy 



luN'erled Hier;u"cliy 



llL,- 



Z5xI0--'eV- 



7.9xlO-5eV- 



i 



7.9xlO'eV 



2.5xJC-%V= 



Figure 1: The relation of neutrino masses and flavor components of mass eigen- 
states. There are two possibilities in the mass relations according to the sign 
of Am23- The case of Atojs > {Arn^s < 0) is called 'normal hierarchy' ('in- 



verted hierarchy' 
relatively large. 



The z/e component of vs is less than 4%. Other mixings are 



As for the Si, the e"^^' term is multiplied by sin 6*13 in the Eq.l^J. Thus 
the effect of CP violation phenomena is proportional to sin^ia. For example, 
the detectability of finite Si from asymmetry of Ve) and P{i'^ v^) 

in near-future accelerator based neutrino experiments depends on the size of 
013 as shown in Fig. |21 Therefore it is important to know the magnitude of 
sin 013 to evaluate the effect of leptonic CP violation, which may explain the 
matter-antimatter asymmetry of the current universe. 

KamLAND experiment measured a deficit of reactor neutrino at average 
baseline 180km. For three generations, what KamLAND observed is expressed 
as, 

P{v^ -> j>e) = cos" 6*13(1 - sin^ 2Bi2 sin^ ^'12) (6) 

There is ~ 5% of ambiguity introduced by the ^13 term, when determining 
sin^ 6*12 from KamLAND oscillation. Therefore precise measurement 6*13 helps 
to obtain more accurate 9i2 from the KamLAND data. 

Thus there are number of reasons why measurement of 0i3 is one of the most 
urgent topics of the neutrino physics. 



1.2 Reactor Neutrino Oscillation 

In reactors, huge amount of low energy D^s are produced in the /3-decays of the 
fission product nuclei. A typical reactor with thermal power of 3GW isotrop- 
ically produces the Avogadro number of j>e's in 15 minutes. The average 
energy detected through i/g + p — > -I- n interaction is around 4 MeV. These 
properties of huge number and low energy of reactor Ve are suitable to use to 
measure neutrino oscillation. 
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Figure 2: Expected 3a discovery regions of sin (5/ as a function of sin^ 20i3 after 
2(f^) and 6.3 (Pfi,) year of exposure in phase II of the JHF-Kamioka project (J- 
PARC). The (black) soHd curve is the case of no background and only statistical 
error of the signal, (red) dashed one is 2% error for the background subtraction, 
and (blue) dotted, (magenta) dot-dashed and (red) solid curves are the cases for 
systematic errors of both background subtraction and signal detection of 2%, 
5% and 10%, respectively fTC| . 
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In 3 generations, the i^e — > Vf. oscillation formula is calculated from the 
Eq.Q and igj as , 

P(Pe = l-4c?3(c?3S?2C?2sin^*2i + s?3C?2sin^*3i+s?3S?2sin2*32)- (7) 

Because |$3i| — |<I'i2 + <&23| ~ |'I'23|, there are essentially two types of oscil- 
lations; one is driven by Am^3 at around 2km and the other one is driven by 
Amj2 around 50km. 

Fig-Elshows the survival probability of 4 MeV in case sin^ 2^13 — 0.1. 




Figure 3: Vg oscillation with 3 generations. 



KamLAND, which measured sin^ 26*12 and Ami2 at around L ~ ISOfcm, 
has demonstrated power of reactor measurement of neutrino oscillation. Fig. 0] 
shows the oscillation pattern of the reactor neutrino measured by the Kam- 
LAND experiment. KamLAND detected an energy-dependent deficit of the 
reactor neutrino. A clear oscillatory pattern is observed in Fig. ^ 

Fig. shows the oscillation parameters measured by KamLAND only (left 
plot) and KamLAND-l- Solar analysis (right plot). From these measurement, 
tan^ e ^ 0.4, Am^ 8 x lO^^eF^ are obtained. 

KASKA aims to measure the sin^ 2^13 at around L ~1.6 km; 1/lOOth base- 
line of KamLAND. At the distance, the Eq.(|2j) reduces to, 

P(Pe -> Pe) = 1 - sin2 2013 siu^ $13 + ©(IQ-^) (8) 

Therefore almost pure sin^ 2^13 can be measured. The current best limit on 
sin^ 26*13 is obtained by CHOOZ reactor experiment in France. The CHOOZ 
experiment used two powerful reactors whose total thermal energy is 8.5 GWt?i. 
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20 30 40 50 60 70 



L„/E_ (km/MeV) 

Figure 4: Oscillation pattern of the reactor neutrinos measured by KamLAND 
at an average baseline 180km. jS] 





Figure 5: Oscillation parameters of electron type neutrinos. The KamLAND 
results (Left) and after combination with solar neutrino experiments. [5j 



7 




Figure 6: CHOOZ sin^ 26I13 limits 



The detector is placed at 1km from the reactors and the result of the measure- 
ments are shown in Fig. They did not confirm the positive result and showed 
the upper limit of sin^ 26*13 < 0.15 for Ato?3 = 2.5 x IQ-^eV"^. 

The Palo Verde experiment also searched for the deficit of reactor neutrino 
with average baseline 0.8km using Palo Verde reactor complex in Arizona, USA, 
whose total thermal energy is 11.63 GW th- They did not observe the deficit of 
the reactor neutrino and the limit of the oscillation measurement is shown in 
Fig.CI 

In view of the currently known Am^^, the baselines of those experiments 
were somewhat shorter than most probable oscillation maximum of 1.8 km. The 
measurement errors were rather large due to uncertainties of reactor neutrino 
flux as well as detector associated absolute errors and low statistics. 

The KASKA experiment is going to measure this oscillation using the D^s 
from world most powerful Kashiwazaki-Kariwa nuclear power station with op- 
timized baseline. Four identical detectors, of which two are placed near the 
reactors and other two are placed at around 1.6 km from the reactors, are used. 
By comparing the data taken at near and far distances, most of the system- 
atic uncertainties are cancelled out and precise measurement becomes possible. 
This far-near cancellation scheme for the ^13 measurement was first proposed 
by Krasnoyarsk group jT2]. Fig. 1221 shows the expected sensitivity for a^-d = 1% 
and 0.5% with 3 years of data. It is possible to obtain 10 times better sensitivity 
than current CHOOZ limit. 
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Figure 7: Result of Palo Verde experiment. JT] 



1.3 Why reactor measurement? 



Near-future long baseline accelerator experiments (LBL) also have sensitivity to 
sin^ 26*13 |15[I13| . In such experiments, high energy f^'s are produced by accel- 
erator and appearance of ly^s is searched for after traveling a few hundred km. 
The probability of j/^ — )■ Ve transition at the oscillation maximum is expressed 
by the following equation. 



f sin2 2013 sin2 0^3 



'12 



Amis 
Am?2 



fg cos 6*13 sin 26*12 sin 2023 sin 20i3 sin 6i 



' cos 9 23 sin 29 12 ^ 



(9) 



where 
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A 



cos(Ai/2) 
1 T AL/iT 
sin(yli/2) 
AL/tt 

V2GFNe, 



and AL/tt ~ 5 x 10^^ for L—300 km. The sign in the definition of / is — 
for the normal (inverted) hierarchy. This measurement depends on many 
parameters which are not necessarily well known. Since sin^ 20i2 and Ar 



"■12 



have turned out not so small by KamLAND and solar neutrino experiments, 
the 2nd term in Eq.@ cannot be ignored any more (this very reason enables 
the measurement of Si in future LBL experiments). Because sin 5; is totally 
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unknown, the 2nd term becomes ambiguity when determining sin^ 2^13 from the 
appearance probabihty. Furthermore, there is an ambiguity of ^23 degeneracy, 
which means there are two possible sin^ ^23 corresponding to measured sin^ 2^23 
value, if it is not unity as, 

. 2„ 1 + Vl-sin2 2023 1- Vl-sin2 2023 
sm^ 6*23 = ^ or ^ (10) 

This degeneracy of ^23 disturbs the ^13 determination by LBL experiments 
further. These situations for the case of sin^ 26*23 = 0.95 and 0.98 are depicted 
in Fig.|Hland|i 




The area bounded by each parabola implies the sinJj-term ambiguity and 
the 023 degeneracy corresponds to two parabolas. Thus even if LBL experi- 
ment measures P(Vfi — *■ v^) precisely, there is ambiguity in determining the 
sin^ 2013. On the other hand, the reactor based oscillation measurement is a 
pure sin^ 20i3 measurement as shown in Eq.®. The expected error on sin^ 20i3 
is ±0.01 and for some parameter regions, the error of reactor measurement is 
better than the ambiguities in sin^ 20i3 determination by LBL experiments. 
Basically LBL and reactor experiment measure different parameters and if both 
results are combined, there are possibilities to solve the 023 degeneracy as shown 
in the Fig. 1101 The possibility to remove the 023 ambiguity by a reactor exper- 
iment was pointed out first in 17 , and was elaborated in detail by KASKA 
group members '18' using a realistic scenario. Thus it is important to know the 
magnitude of sin0i3 to evaluate the effect of leptonic CP violation, which may 
explain the matter- antimatter asymmetry of the current universe. 
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Figure 10: 023 degeneracy can be solved by combining reactor and accelerator 
experiment. 
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Figure 11: Antineutrino running vs reactor-II ((5sin^20i3 ^ 0.01). The 90% 
CL (solid curve) and 3a (dashed curves) aUowed regions in the sin^ 2^13 — 
5;— plane for the true values sin^ 2^13 = 0.1 and Si = 7r/2, are shown. The 
blue/dark curves refer tot the allowed regions for the normal mass hierarchy, 
whereas the red/light curves refer to the sgn(Af3)-degenerate solution (inverted 
hierarchy), where the projections of the minima onto the sin^ 2^13 — (5;— plane are 
shown as diamonds (normal hierarchy) and dots (inverted hierarchy). For the 
latter, the Ax^-value with respect to the best-fit point is also given. Adopted 
from [El. 



12 



KASKA near detectors will collect 600,000 reactor Vf, events in three years 
operation. This will become the highest statistics so far obtained and the ob- 
served reactor neutrino spectrum will become the world standard. The measured 
spectrum can be used to study fission product. Especially there is no calculation 
of reactor neutrino spectrum at energy >8 MeV which are produced generally 
by short hfetime /3-decay isotopes. 

In the future, the reactors v^s may be used to precisely measure 9i2 by 
observing the deficit at baseline 50 km, and Amfg by observing spectral shape 
distortion at baseline ~5 km. Once such parameters are measured precisely, the 
electron type neutrinos can be used as a useful tool to probe internal structure of 
the sun, supernova explosion mechanisms and radiogenic structure of the earth. 
Because Kashiwazaki-Kariwa nuclear power station is the most intense neutrino 
source in the world, it is appropriate to use it to measure such parameters. In 
such an occasion, the KASKA-6'13 detectors can be used as precise near detectors 
as described in the next section. 

2 General Description of the Experiment 

2.1 Kashiwazaki-Kariwa Nuclear Power Station 

KASKA is an abbreviation of ifashiwazaki-ifariwa Nuclear Power Station, fol- 
lowing the tradition to name the experiment after the power station to use. The 
Kashiwazaki-Kariwa nuclear power station 20 , owned and operated by Tokyo 
Electric Power Company (TEPCO), is located 220 km northwest of Tokyo and 
60 km south-west of Niigata-city (Fig. 112(1 . 

The name of the plant comes from the names of the local area on which 
the plant site is locate, Kashiwazaki-city and Kariwa- village. Fig. E| shows a 
picture of the power station. The plant site has an area of 4.2 square kilometers, 
approximately 1.5 km wide and 3 km long. The Kashiwazaki-Kariwa power 
station has 7 reactors which have a total capacity of world's largest thermal 
energy of 24.3GW. The reactors are numbered, from south-west to north-east, 
as #1,^2, #3, #4 and #7, #6, #5. These reactors are arranged approximately 
in a line forming two clusters, each consisting of 4 and 3 reactors. The distance 
between the two clusters is ^1.5 km. Fig.^jcompares the thermal powers of the 
proposed reactor ^13 projects in the world. Kashiwazaki-Kariwa nuclear power 
station is the most powerful neutrino sources, producing an Avogadro number 
of Pg's in every 2 minutes when all the reactors are operating. The reactor #1 
through #5 is the conventional boiling light water reactors (BWR) and #6, #7 
are advanced BWRs (ABWR). Tablenand[21summarize the reactor parameters. 

The newest reactor will be more than 10 years old when KASKA is expected 
to start data taking. Therefore all the reactor fuel is well in equilibrium state 
and average neutrino flux per unit power and energy spectrum are stable. 

The reactor operator TEPCO is very much supportive for basic research. 
KASKA group performed a boring study in the site of Kashiwazaki-Kariwa 
nuclear power station at the site of near detector in the fall of 2004. TEPCO 
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Figure 12: Location of the Kashiwazaki Kariwa nuclear power station. [20] 




Figure 13: Aerial view of the Kashiwazaki Kariwa Nuclear Power Station. [22 
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Figure 14: Site powers for reactor ^13 experiments. 



Table 1: Kashiwazaki-Kariwa Reactor Parameters'^ 



# 


Type 


Thermal Power 
(GWth) 


Start Operation 


1 


BWR5 


3.293 


1985.09.18 


2 


BWR5 


3.293 


1990.09.28 


3 


BWR5 


3.293 


1993.08.11 


4 


BWR5 


3.293 


1994.08.11 


5 


BWR5 


3.293 


1990.04.10 


6 


ABWR 


3.926 


1996.11.07 


7 


ABWR 


3.926 


1997.07.02 



Table 2: BWR, ABWR parameters'^ 



# 


BWR5 


ABWR 


Thermal Power 


3.293(GWth) 


3.926(GWth) 


Electric Power 


1.067(GWe) 


1.315(GWe) 


Fuel Inventry 


132tHM 


150tHM 


Active Core Height 


3.71m 


3.71m 


Active Core Diameter 


4.75m 


5.16m 


Fuel 


UO2 pellet 


UO2 pellet 


Cycle length 


ISmonths 


13months 


Fraction of Core Reloaded Each Cycle 


25% 


25% 


Enrichment 


3.4% 


3.5% 
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and the power station helped much of our activities in the site. 



2.2 Reactor Neutrinos 

In reactors, the power is generated in fission reactions of fissile elements such as 
uranium and plutonium. An example of fission reaction is shown below. 

235u + n ->i40 Te +9^ Rb + 2n. 

The fission products are generally neutron rich nuclei and perform several (3 
decays before they become stable nuclei. 6 /3-decays take place in average per 
fission, which means 6 i^e's are produced per fission along with 200 MeV of energy 
release, resulting in production of 5 x 10^^i7^ per second for Kashiwazaki-Kariwa 
nuclear power station. 

The component of reactor fuel changes along with 'burn-up' through neutron 
absorption and beta decay processes as well as the fission reactions. Below 
shows an example of plutonium production scheme. 

239u ^239 ^ ^ 
239Np^239 pu + e- + ^- 
239pu + n^240 
2«Pu + , 1^241 

^■^^Pu and ^^^Pu produced are also fissile elements and contribute to the power 
and neutrino generation. 238u performs fission reaction after absorbing ener- 
getic neutrons. Fig. 1151 shows an example of change of the fuel components 
caused by the burn-up effect for the case of constant neutron flux. ^'^^U, ^^^U, 
^•^^Pu and ^''^Pu contribute more than 99.9% of the fission reactions. Fig. 1161 
shows the neutrino spectrum for the four main components |23| . Because the 
neutrino spectra depends on fission element, the energy spectrum of the neu- 
trinos from reactor changes along with the burn-up. The fuel exchange takes 
place approximately once per year, in which 1/4 of old fuel rods are replaced to 
new ones. The fuels of all the Kashiwazaki-Kariwa reactors are in equilibrium 
stage, that is 1/4 of the fuel is new and 1/4 of the fuel is 1 year old, 1/4 of the 
fuel is 2 years old and the rest is 3 years old . Thus the component ratio of the 
fuel and neutrino spectrum is stabilized in time. 

In liquid scintillator, i7g's are detected through inverse /3 decay reaction with 
protons. 



i^e + P + n 

The energy threshold of this reaction is 1.8 MeV and the energies of the i7e's 
produced in the /?-decay of U, Pu, Np and neutron produced in the core are 
too low to perform the reaction. Because a typical kinetic energy of recoiled 
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Time evolution of Fissile Elements 




Figure 15: Reactor fuel burn up. 
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neutron is only a few tens of KeV, the kinetic energy of the positron is i^e 
energy minus 1.8 MeV. The positron annihilates with electron and produces 
two 0.511 MeV 7's. Although the lifetime of ortho positronium is 140 ns, the 
pick-off annihilation takes place within a few nano second [^J. Finally, the 
visible energy of the event becomes, 



El^, ^E^- 0.78MeV (11) 

with a good accuracy. Therefore the neutrino energy and visible energy can 
be directly connected and the modulation of neutrino energy spectrum due to 
neutrino oscillation can directly be measured, unlike ve elastic scattering case. 
Because the reaction is the inverse process of the neutron beta decay, the cross 
section is closely related to the neutron lifetime r„, 

27r2SePe ~ , x lO^^^^e (MeV)pe (McV/c) [cm^] (12) 



1.7152m5r„ 



where 1.7152 comes from the phase space factor, including Coulomb, weak mag- 
netism, recoil and outer radiative corrections. Since the of the reaction is 
very small compared with nucleon mass scale, the correction is also small and 
it is known to the precision of 0.2% |25| . 

The visible energy spectrum in case no neutrino oscillation is calculated as, 

/i(i?, t) = Np^a.,{E) J2 ME)n\t) (13) 

i 

where P{t) is the power generation of the reactor at time t, L is the distance 
between the detector and the reactor, a^p{E) is the inverse beta decay cross 
section, fi{E) is the neutrino spectrum per fission of element i, n* is the fission 
rate per unit power generation, and Np is the number of target protons. Due 
to the burn-up effect is a function of time. Fig. El shows the energy 
spectrum of the P^P signals from the four main elements calculated using data 
from ref ■ A typical energy of the detected reactor neutrino is around 4 MeV 
and typical visible energy is 3 MeV. The ratio of event rate is ^^^U: ^^^U: 
^'^^Pu:^'*^Pu~l:l. 39:0. 62:0. 88. Therefore, when calculating expected number of 
events, it is necessary to know the fission rate of each fissile element in the 
reactor cores. Typically the uncertainty of the expected neutrino event rate 
due to the uncertainty of the fuel components is less than 1%. The neutrino 
spectrum from each fissile elements are known in such a way that the neutrino 
detection rate agrees with data at 2.5% of accuracy [211. The thermal power 
generation will be provided from the rector operating company. The error of the 
thermal power measurement is mainly caused by the flow meter of the cooling 
water, which is 2% or less. In total, the absolute error of the expected neutrino 
event is less than 3.4%. However, most of the errors are common for near and 
far detectors and will cancel out when both data are compared and relative 
error on oscillation measurement caused from reactor neutrino flux uncertainty 
is expected to be 0.2% as described in Section |51 
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Figure 17: Energy spectrum for detected neutrino {f^cr,yp) 
parametrization of reactor neutrinos is used I^HI ■ 



A simple 



2.3 General Description of Neutrino Detector 

In KASKA experiment, four neutrino detectors with an identical structure will 
be used. Two near detectors are placed in shaft halls at approximate distances 
of 350'--^400 m from respective reactor cluster and two far detectors are placed in 
a shaft hole at a distance of ~ 1.6 km from all the reactors. Fig.llSlshows relative 
reactor and detector positions. By comparing the data of near and far detectors, 
most of the systematic errors such as detection efficiency and neutrino flux will 
cancel out as will be described in Section|Hl TableOlshows event and cosmic-ray 
rate parameters for KASKA and CHOOZ detectors. Although the KASKA-far 



Table 3: Baseline, event rate and cosmic-ray background ratio. The reactor 
operation, event selection and live-time efficiencies reduce neutrino signal rate 
roughly half. 





target 


Average 


Thermal 


Event 






Detector 


mass 


Baseline 


Power 


rate 


Depth 


ratio 




(ton) 


(km) 


(GW th) 


(/day) 


(mwe) 




KASKA-A 


6 


0.4 


13.2 


520 


90 


1.4 


-B 


6 


0.35 


11.1 


570 


90 


1.6 


-C 


6x2 


1.6 


24.3 


60x2 


260 


1 


CHOOZ 


5 


1.07 


8.5 


40 


300 


0.9 



baseline is significantly longer than CHOOZ one and the depth is shallower, 
the ratio of neutrino flux and cosmic-ray rate becomes better thanks to the 
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Figure 18: Reactor and detector relative positions. 



high power of the power station. Moreover, KASKA detectors have thicker 
buffer layers and the background is expected to be significantly improves from 
CHOOZ. 

Fig. [T^ shows a schematic view of the detector. The central part (Region-I) 
is 6 tons of Gadolinium loaded liquid scintillator (Gd-LS). A candidate LS is 
a Palo Verde type |2Z|, which was successfully used in acrylic containers. The 
liquid scintillator is pseudocumene and mineral oil base with 0.1 % Gd concen- 
tration, with specific gravity of 0.8 g/cm^ and H/C ratio of 1.8. Fig. EHl shows 
a schematic diagram of the neutrino reaction. The reactor is detected in the 
Gd-LS with the inverse (3-deca.y process. 

De + p ^ e+ + 72 (14) 

The positron has kinetic energy of — 1.8 MeV and eventually annihilates with 
an electron, producing two 0.5 MeV 7's as described in the previous subsection. 
The neutron produced has energies of 0(10 KeV) and is thermalized quickly by 
colliding with protons and absorbed by Gd. The Gd produces 8 MeV cascade 
7-rays after the neutron absorption. 

n + Gd^Gd' + 7s(I]£; = 8Mey) (15) 

Because energies of /? and 7 rays of natural radio isotopes are lower than 5 MeV, 
the neutron signal is free from such backgrounds. The neutron absorption signal 
occurs typically 30 after the positron signal. The event is identified by 
requiring that both signals occur within 200 ^s. With this delayed coincidence 
technique, the backgrounds are severely suppressed. 
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Figure 19: Schematic of the KASKA detector. 




Figure 20: Schematic of neutrino signal 
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The Gd-LS is contained in an UV-transparent acrylic vessel with a shape 
of two hemispheres sandwiching a cylinder. The radius of the hemisphere and 
cylinder is 0.9 m and the height of the cylinder is 1.8 m. 

There is a 70 cm thick unloaded liquid scintillator layer (Region-II) outside of 
the Region-I. This region is used to catch the 7-rays which escape from Region- 
I, and to reconstruct the original energies of Gd and positron signals. Fig. |^ 
shows the energy spectra of the Gd and positron signals measured by CHOOZ 
experiment which also has 70 cm 7-catcher region. A clear peak is observed 




Figure 21: Energy Spectrum of (a) Gd signal and (b) positron signal |10j . 

in the Gd signal. Because the event rate at the energy threshold (=5 MeV) is 
small the dependence of the efficiencies on the cut parameter is small. This LS 
will be formulated by mixing PC, paraffin oils and flours, and light output and 
specific gravity are adjusted to same as those of Gd-LS. The region-II scintilla- 
tor is contained in a UV transparent acrylic vessel with the hemisphere-cylinder 
structure whose radius and cylinder height are 1.6 m and 1.8 m, respectively. 
There is a buffer oil layer outside region-II which shields liquid scintillators from 
7-rays and neutrons coming from outside. A PMT support structure divides 
the buffer oil at 90 cm distance from the region-II acrylic vessel. The inner part 
(Region-Ill) mainly works as 7-ray shield coming from the PMT. The outer 
buffer part (Region-IV) is used as a 7-ray and neutron shield from detector 
walls and outer soils, and works as a cosmic-ray veto. In order to improve the 
cosmic-ray veto efficiency, the buffer oil may be made of a weak scintillator 
having 10% of light output of Gd-LS. The base oil will be a paraffin oil having 
a similar specific gravity as region-II scintillator. These oils are contained in a 
stainless steel tank with a diameter 6.5m and shell thickness 15 mm. 
There are 300 10-inch low background PMTs looking at the scintillator region, 
whose photo cathode covers 10% of the surface area. Some LEDs and CCD 
cameras are also arranged in the tank which are used to monitor the position 
of calibration devices. 

Outside of the stainless steel tank is pure water. Water proof PMTs are sub- 
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merged in the water and work as Cherenkov cosmic-ray veto counter. The water 
is contained in a hermetic cyhndrical iron tank. This tank also works as shield 
for geo-magnetic field. Outside of these tanks, there are thick iron layers to 
shield gamma rays from the soils outside. This structure has a leak-safe effect. 
That is, if a crack should develop in the stainless steel tank, the water outside 
would come in the oil region and the oil would not come out. 

There are cosmic-ray tracking devices at the top and bottom of the detector. 
The measured cosmic-ray track is used to estimate spallation background such 
as ^Li and ®He based on the relation of the position of the events and the track. 

In order to reduce the cosmic-ray background, the detectors are placed in 
deep shaft holes with a diameter of 5 m. The depths are at least 150 m for far 
detector 50 m for two near detectors. 

2.4 Sensitivities of KASKA 

2.4.1 sin^ 2^13 measurement 

The expected size of the relative systematic error a^ci on the near-far comparison 
of KASKA experiment is better than 1%, possibly be 0.5% level. Since the 
far detectors will collect 70K reactor j/g with 3 years of data taking and near 
detectors will collect 10 times more data, the statistic error becomes 0.4%. 

The expected sensitivities for rate only and rate + shape analysis are shown 
in the Fig. |^ If no deficit is observed, it is possible to set upper limit on 
sin^ 26*13 at 0.015 to 0.020. This sensitivity is 10 times better than the CHOOZ 
limit. The sensitivity at Am^ > 0.0 leV'^ comes from test of the deficit of De 
flux measured by the near detectors. Because the absolute systematic error and 
statistics error are smaller than previous experiments a sterile neutrino search 
can be performed in this region. 

If the true value of sin^ 26113 is 0.04 or 0.1, for example, and future long 
baseline experiments show Amgi = 2.5 x 10~"^eV^, rate-only analysis constrains 
sin^ 26*13 with 1 degree of freedom (d.o.f.) analysis as shown in Fig. |231 Results 
with rate-|-shape analysis are shown in Fig. 12.4.11 Note that the width of the 
region has the meaning on the line of Am|i — 2.5 x 10~'^eV^ only because of 1 
d.o.f. analysis. 

If a finite sin^ 29i3 is observed by KASKA experiment it means that the 
parameter is reasonably large and there is a good chance to observe leptonic 
CP violation parameter Si from asymmetry of the probability of — > Ve and 

—^ i?e appearances in future accelerator based long baseline experiment. If 
sin^ 26*13 is too small to be detected by KASKA, it means it is difficult to observe 
Si and more sensitive experiments to pin down sin^ 20i3 will be necessary before 
proceeding to accelerator based Si experiments. 

2.4.2 High statistics reactor measurement 

The KASKA near detectors will collect 600,000 reactor events. This will be 
the highest statistics of reactor De detection. Table 0] compares the statistics 
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Figure 22: 90% CL sensitivity to sin'^ 29i3 in KASKA is shown with the analysis 
for 1 degree of freedom. Red curves are of the rate analysis, and thin and bold 
ones are obtained with 0.5% and 1% relative errors in the near-far comparison, 
respectively. 2.5% error on the flux is included in the analysis also. Green curves 
are of the rate+shape analysis with 14 bins within l.SMeV < < 8.8MeV. 
0.5% (1%) uncorrelated relative error between detectors is assumed to be divided 
into 0.5% (1%) uncorrelated error between bins and 0.5% (1%) correlated one. 
The similar assumption is used for 2.5% error on the flux also. The blue curve 
shows the bound by the CHOOZ experiment with the analysis for 2 degrees 
of freedom. The shaded region corresponds to 90%CL allowed region by the 
atmospheric neutrino measurement in SK-I: 1.5 x 10^'^eV^ < Am^ < 3.4 x 
lO-^eV^ 
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Figure 23: Showns are the constraints on sin^ 29i3 to be obtained by the rate 
analysis in the KASKA experiment. The true value of sin^ 2^13 is assumed to 
be 0.04 (red) or 0.1 (green). Amg^ is assumed to be known precisely by the 
future long baseline experiments as Am^i = 2.5 x 10"'^ eV^. The thin and bold 
curves are obtained by (irci — 0.05% and 0.1%, respectively. 
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Figure 24: Showns are the constraints on sin 29i3 to be obtained by the 
rate+shape analysis in the KASKA experiment. The true value of sin^ 2^?i3 
is assumed to be 0.04 (red) or 0.1 (blue). Am|]^ is assumed to be known pre- 
cisely by the future long baseline experiments as Arn^i — 2.5 x lO-^eV^. The 
thin and bold curves are obtained by the different sizes of systematic errors 
which are the same as the values used in Fig. |2S1 
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and detector related absolute uncertainties of various experiments. The abso- 
lute size of the detector related systematic error in KASKA is less than 1.5 % 
thanks to the no fiducial volume cut and low backgrounds. Therefore the data 
taken by KASKA will become standard neutrino Vf, spectrum and contribute to 
other reactor v,, experiments. 



Experiments 


# events (N) 


l/v/iV (%) 


Detector related 
Syst. error (%) 


Distance(m) 


KASKA Nears 


600,000 


0.13 


< 1.5 


350, 400 


cHooz dm 


2,700 


1.9 


1.8 


1050 


Palo Verde HI] 


1,700 


2.4 


4.9 


750, 890 


Bugey EH 


150,000 


0.26 


4.0 


15, 40, 95 


Goesgen [22] 


30,000 


0.57 


6 


37.9, 45.9, 64.7 



Table 4: Comparison of high statistics measurements of the reactor i/g flux. 



A special situation of KASKA is that because 25 % of KamLAND come 
from Kashiwazaki-Kariwa nuclear power station, the KASKA near detectors 
work as near detectors for KamLAND and neutrino flux related error of Kam- 
LAND can partially be cancelled. This effect improves sin^ 20i2 error, together 
with removal of uncertainty of the cos* ^13 term by sin^ 2^13 measurement. 

2.5 Possible Future Extensions of the KASKA project 

Because Kashiwazaki-Kariwa nuclear power station generates the highest D,, 
luminosity in the world, there are a number of opportunities to perform neutrino 
studies with the best precision. 

2.5.1 KASKA-II: Detecting non-0 5i 

In order to proceed to an experiment with errors better than 0.5% level, 12 
tons of the far detector mass is not enough. A 50 tons of detectors and deeper 
shaft holes are under consideration for a next generation KASKA experiment 
(KASKA-II). By combining KASKA-II and data from Super- JHF or Hyper 
Kamiokande, there is a chance to identify non-0 sin J; before proceeding to 
beam if |sini5;| is reasonably large. Fig. 1251 shows the regions in which sin 5; 
is consistent with 0. Therefore, it is possible to identify non-0 sin 5/ if the 
nature chooses the values outside of these regions such as sin^ 26'i3 > 0.05 and 
\5i\ > O.Stt [Sj. 

2.5.2 KASKA-6'12: Accurate sin^ 2^12 measurement 

sin^ 2012 governs oscillation of the electron type neutrinos at large L/E range. 
Natural neutrino sources, such as the sun, the earth, supernovae, etc. are sources 
of electron type neutrinos. Although sin^ 26*12 has been measured by KamLAND 
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Figure 25: sin 5; sensitivity by the combination of KASKA-II, Super- JHF 
(4MW), and Hyper Kamiokande (540kton fiducial mass) |34j. The curves show 
the regions in which data mimics Si = 0. If true values of ^13 and 5i exist outside 
of all of these regions, we can find that Si has non-zero value. The sign of Amgj^ 
is assumed to have been determined by then. 



and solar neutrino experiments, further precision measurement is important 
because it can be used as an inevitable tool to investigate inside of the sun etc. 

Another importance of the precise measurement of sin^ 29i2 appears in the 
relation between the absolute mass of neutrinos and the neutrino-less double 
beta decay searches; The decay is possible for the Majorana neutrino but im- 
possible for the Dirac neutrino. The absolute mass can not be determined by 
oscillation experiments because oscillations depend on Amf^ = mf — rrij but are 
independent of the absolute mass. The information on the absolute mass can be 
obtained by the Cosmological Microwave Background measurements ^|, the 
tritium beta decay measurements and the neutrino- less double beta decay 
searches j32) . The upper bound on the effective neutrino mass of the neutrino- 
less double beta decay can be translated into the upper bound on the absolute 
mass of Majorana neutrinos. The translation, however, depends on the value of 
9i2 very much as we see in Fig. Thus, it is important to know the precise 
value of sin^ 29i2 in order to extract better information on the absolute mass 
from the result of experiments of the neutrino-less double beta decay. If we 
put KamLAND size detector (1 kton) at a distance ~50 km from Kashiwazaki- 
Kariwa nuclear power station, which corresponds to the oscillation maximum of 
Ato^2j it is possible to measure sin^ 26*12 with the highest precision compared 
with solar-pp neutrino experiments. The KASKA-6'13 detectors can be used as 
the precise near detectors. Table El shows a table taken from the reference jSH] 
in which KASKA-6'12 is referred to as the SADO experiment. It shows that with 
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Figure 26: The figure shows the relation between cos20i2 and the upper bound 
on the mass of the hghtest neutrino for a given upper bound (O.leV is used 
as an example) on the effective neutrino mass of the neutrino-less double beta 
decay. The shaded region corresponds the 99.73%CL allowed region of cos 26*12 
obtained by the combined analysis with solar neutrino and KamLAND results 
under the assumption of CPT invariance. 

KASKA-6'i2 configuration, it is possible to measure sin^ 6112 with 2.4% precision. 
In this measurement, background from geo- neutrino is a main issue to obtain 
the highest sensitivity because energy cut introduces large systematics. Note 
that large detector size does not help to reduce geo-neutrino background and 
to obtain a good S/N ratio. One merit of KASKA-6'12 is a good S/N ratio by 
virtue of the high reactor power. Another merit of KASKA-6'12 is that the error 
on the estimation of the geo-neutrino background can be reduced by the help 
of the KamLAND experiment which will have measured the geo-neutrino flux 
of which large part of the sources is common. 



Experiments 


Ss'i^/s'i^ at 68.27% CL 


55^2/^12 at 99.73% CL 


Solar+ KL (present) 


8 % 


26 % 


Solar+ KL (3 yr) 


7% 


20 % 


Solar+ KL (3 yr) + pp (1%) 


4% 


11 % 


54 km 


SADO for 10 GWth-kt-yr 


4.6% (5.0%) 


12.2 % (12.9 %) 


SADO for 20 GWth-kt-yr 


3.4% (3.8%) 


8.8% (9.5%) 


SADO for 60 GWth-kt-yr 


2.1 % (2.4 %) 


5.5% (6.2%) 



Table 5: 
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2.5.3 KASKA-A 771^3: Independent Arri^g measurement 

For the case of three generations neutrino, we have Am^g « ^m^^. However it 
is an important item which has to be tested. If a middle size detector is placed 
at around 5 km from Kashiwazaki-Kariwa nuclear power station, several cycles 
of oscillation can be observed in its energy spectrum. It is possible to measure 

0.1 p . ^ j 
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Figure 27: The solid and dashed lines show the 90%CL allowed regions which can 
be obtained by the KASKA experiment without and with an additional 120 ton 
detector at 5 km from Kashiwazaki-Kariwa nuclear power station, respectively. 
The allowed regions are obtained in the analysis for 2 degrees of freedom in 
order to extract the information on Am§]^. Red and green curves are for the 
case that the true value of sin^ 2^13 is 0.04 and 0.1, respectively. The true 
value of Am|]^ is assumed to be 2.5 x 10~^eV^. The shaded region corresponds 
to 90%CL allowed region by the atmospheric neutrino measurement in SK-I: 
1.5 X 10~^eY^ < Am§2 < 3.4 x lO-^eV^ 

Amjg precisely from the spectrum analysis of the modulation. Fig. shows 
the allowed regions for the positive result in the KASKA experiment with an 
additional 120 ton detector at 5 km from Kashiwazaki-Kariwa nuclear power 
station. The precision of Amf3 is comparable or better than that of the current 
Am23 measurement if sin^ 29i3 is relatively large. Again, the KASKA-6'13 data 
can be used as precise reference spectrum. 

Thus, extensions of KASKA experiment make it possible to measure the 
oscillation parameters of electron type neutrinos by exploiting the large power 
of Kashiwazaki-Kariwa nuclear power station as depicted in Fig. 1281 
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KASKA-High Statistics 
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Figure 28: Scope of KASKA experiments. High statistics measurements of re- 
actor De spectrum with near detectors. High accuracy ^13 measurements at 
L ~ 1.6 km. Precise measurement of Am^g at L ~ 5 km and precise measure- 
ment of 9i2 at L ~ 50 km. 

3 Detector 

3.1 Mechanical Structure 

Fig. El shows a schematic view of the KASKA detector and Fig. 1221 shows 
mechanical structure and size. The detector consists of inner acryhc vessel 
to hold the gadolinium loaded liquid scintillator, outer acrylic vessel to hold 
the 7-catcher scintillator, PMT support frame made of stainless steel, cylindri- 
cal stainless steel tank to hold the buffer oil and cylindrical iron tank to hold 
the water. The shape of acrylic vessels and PMT support frame consists of 
two hemispheres sandwiching a cylinder part. The radius of the hemisphere is 
90 cm, 160 cm and 250 cm for inner acrylic, outer acrylic and PMT support, 
respectively. The height of the cylinder part is common, 180 cm. Fig. 1301 shows 
the structure of the acrylic vessels. The thickness of the acrylic is 1 cm for inner 
vessel and 1.5 cm for outer vessel. The acrylic is a UV transparent type. The 
weight of the acrylic is 0.24 ton for inner and 0.90 ton for outer. The inner 
volume is 7.6 m"^ for inner vessel and 32 m"^ for outer vessel. Compatibility of 
acrylic with liquid scintillator is being tested intensively. The inner acrylic vessel 
is made in a factory and transported to the experimental site and brought down 
to the detector hall. The outer acrylic vessel is divided into two hemispheres 
and two 90 cm- high cylinders. Those units are carried to the experimental site 
by a truck and brought down the shaft hole and assembled in the detector pit. 
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Figure 29: Detector Structure of one unit. 
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Figure 30: Acrylic Vessels 
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There are chimneys with 30 cm diameter (inner vessel) and 60 cm (outer vessel) , 
which are to be used to put calibration devices in and to control the liquid level. 
The chimneys are made of black acrylic. 

The PMT support frame has a jungle-gym type structure. The number of 
PMT is ~300 and 1.5 tons of buoyancy force of the PMTs will pull the structure 
up when the oil is full. The PMT support structure has counter weights to cancel 
this buoyancy force. 

The diameter of the stainless steel cylindrical tank is 6.5 m and the height is 
990 cm and the the volume is 330 m'^ . The thickness of stainless steel cylindrical 
tank is 15 mm and its total weight is 32 tons. This tank holds 260 tons of oils. 
The tank has to strong enough to self-stand when the water around is empty 
while endure the buoyancy force of ~ 20 ton when the water is full. 

The diameter of the cylindrical iron tank is 7.5 m and the height is 1,060 cm 
and the inner volume is 470 m^. The thickness is 11 cm to shield 7-rays coming 
from outside. It is considered to make this iron tank as geo-magnetic shield by 
de-magnetizing the iron. Table El summarizes detector parameter structures. 



Table 6: Detector Structure Parameters. 



Item 


R 

(cm) 


H 

(cm) 


t 

(cm) 


V 

(m3) 


W 

(ton) 


Inner acrylic 


90 


180 


1.0 


0.20 


0.24 


Outer acrylic 


160 


180 


1.5 


0.75 


0.90 


PMT support 


250 


180 






2 


S.S. tank 


325 


990 


1.5 


4.0 


32 


Iron tank 


375 


1060 


11 


37 


290 



R: Radius of cylincr/hcmisphcrc, H: Height of cylinder/cylinder part, t: Thick- 
ness of the vessel, V: Volume of the structure material and W: Weight of the 
vessel. 

3.2 Liquid Scintillators 

The KASKA detector uses four kinds of liquid. From inner to outer, v target 
scintillator (Region-I), 7 catcher scintillator (Rcgion-II), buffer oil (Region-Ill, 
IV) and water for cherenkov cosmic-ray anti-counter (Region- V). Tabled shows 
components of each candidate liquid. 

In the table, BC521 is a commercial liquid scintillator of a Saint-Gobain 
product with Gd concentration of 0.5%. PC stands for pseudocumene (1,2,4 
tri-methyl-benzene (C6H3(CH3)3) ), TD stands for Tetradecane (C14H30). v% 
means % in volume and w% means % in weight. Table |S1 shows properties of 
these liquids along with the mixed ones. 
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Table 7: Liquid Scintillator Components 



Region 


Mixture 


Volume (m'^) 


Weight (t) 


I 


PC(20v%)+TD(60v%)+BC521(20v%) 


7.6 


6.1 


II 


PC(15v%)+TD(85v%) 


24 


19 




+PPO(0.6w%)+bis-MSB(0.01w%) 






III+IV 


Paraffin oil (+PPO) 


300 


240 


V 


pure water 


160 


160 



Table 8: Liquid Properties 



liquid 


mass 


specific 


flash 


H/C 


light 






gravity 


point 




output 




(ton) 


(g/cm^) 


rc) 




(% anthracene) 


BC521 m 




0.89 


44 


1.314 


68 


PC 




0.880 


54 


1.33 




TD 




0.763 


99 


2.14 




I 


6.1 


0.79 


49.5 


1.8 


53 


II 


19 


0.78 


76 


2.0 


53 


III+IV 


240 


0.76 


99 


2.14 


5 


V 


160 


1.00 
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3.2.1 Neutrino Target (Region-I) 



The active neutrino target is Gd Loaded liquid scintillator (Gd-LS) with Gd 
concentration of 0.1 %. A candidate is a Palo Verde type Hquid scintillator |S7). 
This LS is formulated by mixing 20 % of BC521, 20% of PC and 60% of TD. 
The molecular formula of the Gd in BC521 is 'gadolinium 2-ethylhexanoate'; 
Gd(CH3(CH2)3CH(C2H5)C02)3) xHaO]. QJx- The light output is 53 % an- 
thracene. Along with 10 % of PMT cathode coverage with its 20 % quantum 
efficiency, the light yield will be 150 p.e./MeV, which corresponds to statistical 
energy resolution of 

^ . (16) 

E ^S(MeV) ^ ' 

H/C ratio is 1.8 and proton density is 7.7 x 10^^/m^. The Palo Verde experiment 
used this kind of LS contained in 7.4 m long acrylic cells for a few years without 
problems. TD is chosen as the base solvent because it is mono- molecular parafhn 
oil and proton density is exactly known. Paraffin oil does not have double 
bond in its molecular structure and does not absorb scintillation light, and it 
is immune to oxidization. Because TD and PC are mono-molecular material, 
the temperature window of the distillation in refinery process is narrow and 
contamination is intrinsically small. This strategy was very successful in the 
KamLAND liquid scintillator. TD has higher flash point and slightly larger H/C 
ratio and specific gravity than dodecane which is successfully used in KamLAND 
detector. 30 m^ of the Gd-LS for four detectors will be formulated at the 
same time and stored in a storage tank in order to make the components of 
the four detectors same. When filling in the detector, the LS is carried by 
200 liter stainless steel drum cans. The weight of each drum can is measured 
with relative error of <0.3 % using carefully inter-calibrated scales before and 
after emptying into the detector and the total amount of Gd-LS introduced is 
measured precisely. The acrylic vessel containing the Gd-LS has a chimney with 
30 cm diameter. 1 cm of head level difference corresponds to 0.01% of LS volume 
difference. By monitoring the head level, precision monitoring of LS volume in 
the fiducial region can be performed. 



3.2.2 7-ray catcher (Region-II) 

The purpose of this liquid scintillator is to detect the escape gamma rays from 
region-I and reconstruct original energies of positron annihilation and neutron 
capture on Gd signal. The light output is adjusted to be the same as Region-I 
LS. A candidate of this LS is, 

PC(15%)+TD(85%)+PPO(0.6w%)-f^bis-MSB(0.01w%) 

The specific gravity is slightly smaller than Region-I LS to keep safety of the 
acrylic vessel. The concentration of PC is reduced from that of Region-I LS to 
raise the flash point for safety reason. 
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3.2.3 Buffer oil (Region-III&IV) 



The purposes of the buffer oil is to shield the liquid scintillators from 7-rays 
from the PMT glass and outside soils, and to tag the cosmic-rays. 300 m'^ of 
paraffin oil will be used as the buffer oil. In order to detect stopping muons 
whose energy is less than the Cherenkov threshold and to improve cosmic-ray 
veto efficiency a slight light output 10 % of Gd-LS light output) may be 
given to the oil by adding PPO. The specific gravity is slightly smaller than 
Region-II scintillator. 



3.2.4 Cosmic-ray veto (Region- V) 

Outside of the stainless steel spherical tank is filled with 160 m'^ of pure water. 
The purpose of this region is to tag the cosmic-rays to veto the cosmic-ray related 
background and is to shield the 7-rays and neutrons from the soil outside of the 
cave and concrete to form the cave. The water is kept purified by circulating. 



3.3 Photosensor 

In order to detect the scintillation photons from liquid scintillators, large aper- 
ture photomultiplier tubes (PMT's) are used. One of the candidates is Hama- 
matsu Photonics K.K. (HPK) 10 inches PMT (R7081). In order to obtain 
necessary photoelectron yield, more than 10% of the photo cathode coverage is 
necessary and >300 10 inch PMT's are necessary. Fig. -|^ shows the shape of 
the PMT and table IHI shows its specifications. Fig.-|^also shows the quantum 




Figure 31: Shape of the PMT envelop and quantum efficiency. jHl] 

efficiency for the PMT. The photo cathode is bi- Alkali type with 25% quantum 
efficiency around 400 nm wavelength. 

Single photoelectrons response for the PMT was measured as follows. As a light 
source for this measurement, we used a PiLas picosecond laser light system (= 
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Table 9: PMT specifications ^ 



Type No. 


minimum 


Approx. 


Q.E. 


Rise 


TTS 




effective area 


Weight 


at 390nm 


Time 






(mm) 


(g) 


(%) 


(ns) 


(ns) 


R7081 


220 


1150 


25 


4.3 


2.9 



438.7 nm) having a high precision pulse timing (<30 ps). In order to obtain a 
single photoelectron response we used neutral density filters to reduce the laser 
light. For the digitalization of the pulse amplitude we used a LeCroy 2249 A 
ADC. Fig. 1321a) shows the single photoelectron spectrum for R7081. A clean 
single photoelectron peak can be seen at around 70 ADC channels. Fig. 1321b) 
shows the leading edge timing for the single photoelectron signals. We used a 
TDC (Houshin Co. Ltd) for the timing measurement. After a time walk correc- 
tion using the pulse height information (assuming a square root dependence on 
the pulse height), a transit-time -spread (TTS) has been obtained to be around 
2.5 ns, which is better than the catalogue value of 3.7 ns in which the time walk 
correction may not be applied. Fig. |321 c) shows the timing distribution after 
the time walk correction. 

One of our worry for the large size PMT's is response uniformity over the 
large surface area. By changing the laser light incident position, we have mea- 
sured the pulse height spectra for three different positions. Fig. ESI a) shows the 
peak ADC channels for single photoelectron spectra measured at three incident 
positions: (center), 5 cm and 10 cm from the PMT surface center. Data for 
R5912 are also shown for comparison. Both R5912 and R7081 have a good 
uniformity within 10 % variance. Fig. 1331 b) shows the uniformity of the res- 
olutions (standard deviations). They are also very uniform over the measured 
surface area. We can conclude from the above measurements that both R5912 
and R7081 have enough performance for the detection of the scintillation light. 

Generally PMT glass contains significant amounts of radioactive elements 
and causes single background hits. In order to reduce this background, a low 
background glass is developed. Table [W\ compares the U/Th/K contents of the 
PMT glass used in various experiments. With this glass and 90 cm of buffer 
thickness between glass and region-II scintillator, the single rate of the detector 
will be a few Hz at above 0.7 MeV threshold, which is satisfactory value. 

There is an optional idea to install some 2" PMTs besides the 10" PMTs to 
measure the total energy of cosmic ray signals without saturation problem. 

Fig. 151 shows the PMT positions. 300 PMTs are uniformly distributed. 

3.4 Cosmic-ray Tracking Device 

A cosmic-ray tracker system will be installed in the detector. The role of the 
cosmic-ray tracker is to track the cosmic muon passage which goes through 
the detector. Spallation background events are estimated by its association to 
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Figure 32: a). Pulse height spectrum and b) leading edge timing distribution 
for the single photoelectron signals obtained with R7081. c) Timing distribution 
for R7081 after the time walk correction. The TTS of 2.5 ns is obtained from 
this figure. 



Table 10: radio purity of PMT glasses 



Type 


my 


^^^Th 


4UK 


Reference 




(ppm) 


(ppm) 


(ppm) 




KASKA candidate 


0.046 


0.12 


0.01 




CHOOZ 


0.07 


0.1 


0.018 




SuperK/KamLAND 


0.48 


0.47 


0.08 


41, 


Kamiokande 


0.25 


0.32 


0.016 


41 
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Figure 33: a) Uniformity of the response for three R5912 and one R7081. The 
incident positions of the laser hght were at (center), 5 cm and 10 cm from the 
PMT surface center. Peak values were stable within 10%. b) Uniformity of the 
resolution for three R5912 and one R7081. 
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Figure 34: PMT positions. Left: top/bottom view of the hemisphere part. 
Right: unfolding view of cylinder part. 
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the cosmic-ray track. The position resolution required is a few centimeters, 
comparable to or better than the resolution of the vertex reconstruction for the 
neutrino events in the target region. 

To achieve these purposes, several options can be considered for the hardware 
implementation: use of liquid-scintillator(LS)-based, plastic-scintillator(PS)-based 
or gaseous detectors such as RPC. In the case of a gaseous detector, building a 
gas-supplying system down in the shaft hole would be complicated, and some 
safety issues might arise, too. Not much attention has been paid on this option 
at this moment. 

In cases of LS-based or PS-based detectors, the detector will consist of scin- 
tillator strips with a segmentation to match the position resolution. In a design 
investigated, the scintillator is a rectangular with 10 cm in width, 3 cm in height, 
and about 5 m in length. In order to cover the whole area of the main detec- 
tor, about 50 scintillator strips are arranged. A layer consists of two planes, 
X and Y, which give a 2-dimensional position. Fig.-|2Sl shows the scintillator 
layer. These two layers are set up on the top and bottom of the main detector. 
In reality, there have to be some dead spaces to accommodate the mechanical 
support structure. 




Figure 35: Schematic drawing of a layer of cosmic-ray tracker. 

A wavelength-shifting (WLS) fiber is used to transport the scintillation light 
to the photodetectors such as PMTs. By using the WLS fiber, the photodetector 
can be set up in a convenient place. The scintillation light is read from both 
sides. The time difference between both sides will also give a redundant position 
measurement. The WLS fiber considered is about 1.5 mm in diameter, and 
a multi-clad type. Because its angle of the total reflection compared with a 
single-clad type is larger, the light yield is larger as shown in Fig.-ISSI It is 
necessary to consider the attenuation length of the WLS fiber. For example, the 
attenuation length of a multi-clad type fiber, Y11(200)MS made by Kuraray, is 
about 3.6 meters. The hght through the WLS fiber is read by Multi- Anode PMT 
(MAPMT) such as Hamamatsu H7546 series. The advantage to use MAPMTs 
is to save the cost, and a smaller occupied space compared with PMTs installed 
for all cells. 
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Figure 36: Sketch of Y-11 WLS fiber. 



3.4.1 A case of LS-based detector 

The advantage of this option is low cost scintillator. An aluminum or SUS 
case and an acrylic cap can be used for the container. UV-transparent acrylic 
should be used for the cap. There is a ditch on an acrylic cap, and the WLS 
fiber is inserted as shown in Fig.-|23 In this case, it is necessary to consider 
the attenuation length and light yield because the solid angle of the fiber to 
absorb the scintillation light is small. Multiple readout fibers may be necessary 
to collect enough photons. 




Figure 37: Drawing of an acrylic cap and inserted WLS fiber at the end of a 
LS-based detector. 



3.4.2 A case of PS-based detector 

In this case, extruded plastic scintillator ^38 can be used. This scintillator has 
inferior optical property, but costs much less than usual PS. The scintillator 
will be painted on the surface to reflect the light, and a hole will be made to 
pass the WLS fiber. Extruded PS is used by the SciBar detector at K2K02j 
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and by the MINOS experiment [?]. Basic characteristics of this scintillator is 
shown in table ^2 The structure of the plastic scintillator is shown in Fig. |2HI 
WLS fiber is inserted in the center of the scintillator. The advantage is that it 
is not necessary to consider the attenuation length of scintillator. However, it 
is necessary to consider the attenuation length of the WLS fiber compared with 
the case of a LS-based detector. 
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Figure 38: Drawing of an extruded PS strip. The unit is mm. 



Table 11: Basic characteristics of the scintillator used in SciBar detector 



Scintillator material 


polystyrene with PPO (1%) and POPOP (0.03%) 


Emission wavelength 


420 nm (blue) 


Dimensions 


2.5 X 1.3 X 300 cm^ 


Hole diameter 


1.8 mm 


Reflector material 


Ti02 (15%) infused in polystyrene 


Reflector thickness 


0.25 mm 



4 Shaft Hole 

The required accuracy on sin^ 20i3 measurement is a fraction of percent. It is 
necessary to suppress the cosmic-ray related backgrounds significantly to realize 
such accuracy by placing the detectors deep underground. For KASKA case the 
far detector should be placed at least 250 meter water equivalent (m.w.e.) deep 
to sufficiently suppress cosmic ray muons and near detectors should be placed 
at around 80m. w.e depth. Because there is no mountain around Kashiwazaki- 
Kariwa nuclear power station for such distances less than 2 km, KASKA uses 
vertical shaft holes. As the density of the soil in this area is 1.75 g/cm"^, the 
physical depth of the shaft holes are 150m for far detector and 50m for near 
detectors. Fig. 123 shows the schematic of the shaft hole and detector pit under 
consideration. The diameter of the shaft hole is 5m. A 18m long horizontal 
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Figure 39: Dimension of a shaft hole for the far detector (one detector case). 
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corridor connects the shaft hole and the detector pit. The detector pit has a 
cyhndrical hole of 8m diameter and 10m height. The neutrino detector will 
be placed in the pit. The construction method of the holes will be a short- 
step method. A few meters are dug down at one time by blasting, then debris 
are taken out by a crane and then the wall is reinforced by concrete. These 
procedures are repeated until the hole reaches the target depth. Fig. ^ shows 
an example of how such construction proceeds. Constructions of such a vertical 
shaft hole have been performed in several locations. For example, Japan Atomic 
Energy Association has been digging l,000m-depth shaft holes with diameter of 
4.5m and 6m in parallel for the purpose of studying deep underground storage 
of nuclear waste. Fig. 0| shows a figure of the twin shaft holes and a picture 
taken while a shovel car is cleaning the debris. R&D for digging the KASKA 
holes has been under way. Extensive studies of the geology were performed 
by the electric power company when the reactors were constructed and overall 
structure of the geology around this area is well known. Fig. 1421 shows stratum 
at near-B detector position which was measured by our boring study. Table IT^ 
shows the properties of strata. 



Table 12: Properties of strata 



Name 


main constituents 


density 
(g/cm^) 


water component 
(%) 


Shinki-Saso 


fine sand 


1.64 


5.6 


Bansya-Saso 


rough sand 


1.83 


15.1 


Yasuda-so 


silt, clay 


1.76 


47.0 


Nishiyama-so 


silt rock/clay rock 


1.74 


47.4 



The general property of the geology is low density mud with abundant wa- 
ter component. A boring study for the KASKA-B detector was performed in 
fall 2004. Fig. 01 show a picture taken during the boring study. Sample cores 
were taken from several depths and properties of the soils, as well as U/Th/K 
contaminations, have been studied. The cosmic-ray and 7-ray background were 
also measured at various depths using Nal and plastic scintillators. These in- 
formation will be used for detailed design of the KASKA neutrino detector. 

5 Readout Electronics 
5.1 Property of Signals 

The KASKA experiment is composed of four identical detector units each having 
300 photomultiplier tubes (PMT), to view the Region-I and II liquid scintilla- 
tors. Because there are 4 detectors, the total number of PMTs is 1,200. Those 
signals from the PMTs are read out by flash ADCs (FADC). There also will 
be hundreds of signals from cosmic-ray anti-counter and cosmic-ray tracking 
devices. Signals from those devices are readout with conventional ADC/TDC 
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Figure 40: Example of short step excavation procedure. Prom left to 
right, drilling/setting an explosive, blasting, taking out debris, moving 
down/ concreting. 
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Figure 41: 1000 m depth twin holes and its construction sciene |^ 
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Figure 42: Stratum at near-B detector position, measured by boring study. 
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systems. The following description concentrates on the readout electronics for 
the central PMT data, which are being developed for this experiment. Signal 
sources are anti-neutrino interaction, cosmic rays and other backgrounds and 
calibration signals. Detection and identification of anti-neutrino is performed 
by a prompt and a delayed signal generated in an inverse /3-decay as described 
in section The former signal is produced by a positron and the latter is 
created by 7 rays from gadolinium nucleus after absorbing neutron. Fig. 1441 
shows the conceptual picture of the Vf, signal. 
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Figure 44: Characteristics of signal. 



The energy of the prompt signal is 1 to 9 MeV and that of the delayed signal 
is 8 MeV. The mean interval time of the two signals is about 30 /iS. The rates of 
neutrino signals per detector are ~l/hour for far detector and r^lQ/howc for near 
detector. The dominant background source is cosmic-ray muons. The maximum 
energy deposit of the cosmic-ray in the scintillator is 800 MeV for single muon 
and much more for showering muons. The rate of the cosmic-rays are Hz 
for far detectors and ~170 Hz for near detectors, respectively. The cosmic-ray 
signal has to be separated from the neutrino signal with very high efficiency 
for the suppression of backgrounds. An important effect to be considered is 
production of radioactive elements with spallation reaction. It leaves some radio 
activity which remains after the cosmic-ray passing. These activities have to be 
eliminated from the neutrino signals also. As the energy deposit of the cosmic- 
ray signals is huge compared with the reactor neutrino signals , the electronics 
system will be saturated. It is, therefore, necessary to realize fast recovery from 
cosmic ray signals for the deadtime-less measurement. In addition, there are low 
energy single background signals which come from radioactive elements in the 
detector material and soils outside of the tunnel, such U and Th, and their decay 
chains and also from and ^°Co, whose visible energy is less than 5 MeV. It 
is, therefore, expected to be reduced to less than 10 Hz with 0.7 MeV energy 
threshold by putting a thick absorbing buffer. Finally, the positron produced 
by stopping muon (Michel positron), whose maximum energy is 53 MeV, will 
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be used for the calibration in the energy range up to a few 10 MeV. 

/i+ e+ + + j^/i 

As the hfetime of muon is 2.2/^s, the electronics has to recover quickly from 
the saturation caused by the initial signal. 

Under these circumstances, the requirements for electronics are as follows. 

• Measurement of output charge of each PMT hit to reconstruct energies of 
the signal. 

• Measurement of the timing of each PMT hit to reconstruct positions. 

• Sufficiently wide dynamic range for neutrino signals and cosmic ray signals. 

• Dead-time-less Multi-hit readout capability. 

• Measurement of the timing separation between two contiguous events. 

• A fast recovery from the saturation of large cosmic-ray signals and moni- 
toring after- muon activities. 

In KASKA experiment, most of those information are measured by FADCs, 
with multi-bank memories for deadtime-less measurement with precise time 
stamps used to measure the time interval of the prompt and delayed signals. 

5.2 Electronics Scheme 

Fig. 031 shows a schematic diagram of the readout electronics. A signal from 
the PMT is amplified and put into a 500MHz FADC together with the marker 
or calibration signals. The output from the FADC is divided into a storage 
circuit and a trigger circuit. The FADC data going to the storage circuit are 
recorded in a bank memory of 300 ns equivalent depth, after a digital delay 
which compensates the latency of the trigger decision. There will be two, or 
more if necessary, bank memories to eliminate the dead time. When a trigger 
is issued, the direction of the FADC data flow is switched from a bank to the 
next one. The newly coming data being digitized are recorded/shifted in the 
new bank memory, while the stored data in the first bank memory are being 
readout by a computer. Each bank has a register which holds the data of 
absolute timing when the trigger is issued. The FADC data for triggering is 
digitally time-integrated to get the total charge information of the signal. If 
the total charge exceeds 0.25 p.e. equivalent, the PMT is considered to have a 
signal and a PMT trigger is issued. This scheme has been provn to work with 
latency of ~80 ns for a prototype. Then a PMT trigger information is sent to a 
trigger decision module, where the number of hits from all the PMTs is counted 
up. If the number of hit PMTs exceeds 0.7 MeV threshold, a global trigger is 
issued, which is sent to a trigger control circuit. The trigger control combines 
the amplitude and timing information of the data and decides if the currently 
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Figure 45: Schematics of KASKA electronics. 
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active bank is switched and send information about which bank to read out to 
the data acquisition computer. 

Fohowings are assumed KASKA detector parameters and the properties of 
signals for the rough estimate of necessary gain and dynamic range which are 
required for the electronics. 

1. KASKA detector parameters 

• Required energy loss for one photon in the scintillator: 150 eV 
=> 46% Anthoracene equivalent. 

• PMT photocathode diameter: 22 cm. 

• Radius of photocathode position in detector: 3 m 

=J> photocathode coverage=10 %. 

• Radius of the region-I scintillator: 1.4 m. 

• Quantum efficiency of PMT: 20 %. 

• amplification gain of PMT: 5 x 10^. 

• Rise/Fall time of single photoelectron signal: 10/15 ns. 

• Specific impedance of the signal cable: 50 fl. 

• Load impedance of PMT output including a back termination: 25 SI. 

• Samphng rate of the FADC: 500 MHz. 

• Input voltage range of the FADC: ±0.3 V. 

• Resolution of the FADC: 8 bits => 2.3 mV/ch. 

• Amplifier gain: 10. 

• Attenuator ratio: 5:1. 

• Trigger threshold for one PMT: 0.25 p.e. equivalent. 

2. One photoelectron signal 

Approximating the signal shape to a triangle with 10 ns rise time and 15 ns 
fall time, the peak voltage of the one photoelectron signal is 1.6 mV. After 
the amplifier of gain 10, it corresponds to 7 FADC counts. The signal is 
distributed to 13 bucket wide. The trigger discrimination is applied to 
the output data of digital integrator. The integrated FADC counts for 
one photoelectron is 44, and 0.25 p.e. PMT trigger threshold corresponds 
to 11 counts. In actual signals, the shape of the tail is exponential and 
the signal will distribute to much wider time span than 25 ns and the 
integration period will be 50 ns. 

3. Low energy signals 

An average number of photoelectrons created at one PMT for 1 MeV 
signal is 0.45 p.e. Thus half of PMT do not have a hit and suppression 
is effective to reduce the data size. The maximum energy for the reactor 
neutrino signal is 9 MeV, which produces 4.1 p.e. /PMT of scintillation 
light if the event occurs at the center of the detector. Because events 
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occurred near the acrylic wall deposit 3.5 times more light to the nearest 
PMT than those occurred at the center of the detector, the nearest PMT 
will generate 14 p.e./PMT for such peripheral events with energy 9 MeV. 
The dynamic range of the FADC corresponds to 37.5 p.e., which safely 
covers the maximum reactor neutrino energy. 

4. High-energy signals 

A penetrating cosmic muon causes as much as 750 MeV of energy deposit, 
because the maximum track length in the scintillators is 5 m. 750 MeV 
of energy deposit produces 350 p.e./PMT of scintillation light. This is 
out of dynamic range of the FADC system and at least 1/5 of attenuation 
is necessary. If the scintillation light comes at the same time, the pulse 
height from the PMT becomes 1.5 V. The output of the amplifier will 
saturate for this kind of signals and fast recovery from the saturation is 
required to the amplifier. A clipping diode can be installed to the amplifier 
input. 

Another idea to treat this issue is to install smaller diameter PMTs in 
addition to 10" PMT, which is now our baseline scheme. If 2" PMT is 
used, the aperture of the PMT is 1/25 of the 10" PMT and the number of 
photoelectron obtained by PMT is as much smaller and the PMT will not 
saturate even for cosmic-ray signals. This information is used to judge the 
cosmic-ray is showering or not and set additional dead time for showering 
muons. 

5.3 Trigger 

The basic trigger is issued for events which have energy greater than 0.7 MeV. 
All the events with the energy greater than the threshold are recorded. For 
this energy the number of expected photoelectron is 0.3 p.e./PMT. Thus the 
energy of the signal approximately corresponds to the number of hit PMTs and 
the global trigger is issued based on the multiplicity of the hit PMTs. The hit 
multiplicity for 0.7 MeV signal is about 90 PMT hits. A PMT hit is defined 
by 'PMT trigger', which is issued when the 50 ns-wide integrated FADC hits 
becomes 0.25 p.e. level or more. The PMT hit information is also used for 0- 
suppression procedure. The integration of the FADC data and the summation of 
PMT hits are performed by field programmable gate arrays (FPGA) as shown 
in Fig. 221 The information of the PMT trigger from all the PMTs are sent 
to the global trigger module and total number of hits with integration time of 
100 ns is calculated. The trigger circuit can generate more sophisticated trigger, 
such as delayed coincidence -t- low threshold trigger, pre-scaled trigger, random 
trigger, GPS trigger for absolute time stamp, forced trigger for calibrations, etc. 

5.4 Readout modules 

The FADC module is composed of a daughter board, which is a main part with 
a FADC chip, and a mother board which works mainly as an interface to the 
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readout bus. 



1. Daughter board: It includes a buffer amplifier wliose gain is selected either 
10 or 1/5 by a jumper. The input impedance is selectable either high or 
low by a jumper, so that two FADC units can be connected by a daisy 
chain of one input as a high gain FADC and a low gain one. Fig. BHl 
shows a schematic of the input circuit of FADC module. It is also usable 
as independent FADC channels, as well. The amplifier is followed by a 
FADC chip, which is 8-bit 500 MSPS ADC as a baseline. It wiU be driven 
by a 500 MHz clock in the present design. The input voltage range is 



2. Mother board: It includes FPGAs which work as controllers of the module 
and an interface with rear-end electronics which will be either Compact- 
PCI or VME. 



Timing depth of one peace of data is 300 ns, which is determined by both 
pulse duration for the fast component and the expected latency of the trigger. 
The FADC chip is always active and the digitized data are written on a bank 
of memories in a FPGA. When the module receives a trigger from the global 
trigger circuit, the currently active bank and associated counter is frozen and 
the module activate the next bank in turn. The frozen memories are then read 
out to the rear-end electronics, while the next data are independently read into 
another bank. This feature of multi-buffer memories is adopted for deadtime-less 
DAQ, that is essential for the present experiment. The time interval between a 
buffer and a successor is measured by the absolute time stamp. [H] 



±0.32 V. 




Figure 46: FADC modules input 
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6 Calibration 



6.1 Introduction 

Because reactor ^13 measurement needs a high precision and rehabiUty, intense 
calibrations arc necessary throughout the run time. CaUbration sources are clas- 
sified into two categories. One category is optical and electronic sources, such as 
diffusive laser ball, optical fibers, LED and so on. The sources are used to mea- 
sure responses of PMT and electronics and performs relative calibration. The 
other category is energy and position calibration sources. These are radioactive 
sources (gamma source, neutron source, positron source). Gamma sources are 
used for determination of absolute energy scale and energy resolution. Neutron 
and positron sources are used for measurements of detection efficiency and re- 
sponse. In addition, in order to locate a source at aimed positions, a calibration 
manipulator system will be developed. 





Source 


Feature 


Optical 


Diffusive laser ball (SNO) 


Isotropic laser 




UV laser (CTF-Borexino) 


Pulsed laser 




LED 


simple 


Neutron source 


^^^Am-^Be 


with 4.4 MeV gamma ray 






spontaneous fission 


Gamma source 




0.662 McV 




^"Mn 


0.835 McV 




*^"Co 


1.173, 1.333 McV 






1.35 MeV 






1.5 MeV 




H neutron capture 


2.22 MeV 




^"^C neutron capture 


4.95 MeV 




Gd neutron capture 


8 MeV 




p-T source (SNO) 


19.7 MeV(Q value) 


Positron source 


-^Na 


2.8y. 0.511 MoV. 1.275 MoV 




''^Ge 


280d, alternative energies 




Michel positron 


Emax=53 MeV 



Table 13: various sources 



6.2 Optical and electronic calibrations 

The purpose of optical and electronic calibrations is to check the response of 
each PMTs and to monitor the attenuation length of liquid scintillators. An UV 

pulsed laser can be used for monitoring of performances such as one p.e. gain, 
uniformity, timing and pedestal, of each PMT. An isotropic light source like a 
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diffusive laser ball (SNO) can be used for monitoring the attenuation length of 
the detector components. 

6.3 Energy calibration 

In an antineutrino event, the positron produced by inverse-/3-decay interaction 

loses its energy and annihilates with electron in LS and the total energy deposit 
is detected as a prompt signal by the PMTs. In addition, neutron produced by 
the inverse-beta interaction is captured by Gd after about 30/xsec, and gamma 
rays are emitted that are detected as a delayed signal by the PMTs. Therefore 
7 ray source, positron source and neutron source are used for energy calibration. 



6.3.1 7 ray sources 

The purpose of using 7 ray sources for calibration is to determine the absolute 
energy scale and energy resolution. Gamma ray sources can be classified into low 
energy region corresponding to the prompt signal response (0.7^9.0 MeV) and 
high energy region corresponding to the delayed signal response (5.0~11 MeV). 
There are many kinds of gamma ray sources in the low energy region. For 
example, i37Cs(0.662MeV) and 6nCo(1.173McV+1.333MeV) arc useful sources 
in the low energy region. On the other hand, there are little kinds of gamma 
ray sources in the high energy region, because it is difficult to produce long 
life and high rate sources. Therefore, if we need gamma source in high energy 
region, one option is to develop an original source like pT source used by SNO 
experiment. 

6.3.2 Positron sources 

The purpose of calibration using positron sources is measurement of positron 
response. ^^Na is a very useful source because of having longer life (2.8 years) 
than other positron sources. Positron produced by decay annihilates with 
electron in the source package, and then two gamma rays (0.511 MeV) emit 
back to back. In addition, one more gamma ray (1.275 MeV) emits. These 
annihilation gamma rays are used as mimic signals. ^^Ge is a useful source, too, 
because of emitting gamma rays with alternative higher energies, but the half 
life is relatively short (287 days). Another way is to use cosmic rays. If stops 
in the LS, it decays and produce energetic positrons in the LS. 

At+ ^ e+ + i/e + i/p. (17) 

The energy distribution of this positron (Michel positron) is well known and it 
can be used to calibrate energy range up to 50 MeV. 
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6.3.3 Neutron sources 

The purpose of calibration using neutron sources is the measurement of neutron 
capture response by the nuclei (H, ^'^^Gd, ^^''Gd) in LS. Gamma rays produced 
by neutron capture of H, ^^^Gd and ^^^Gd give a total energy of 2.223 MeV, 
8.536 MeV and 7.937 MeV respectively. Especially, the neutron capture by 
""^^^Gd and ^^^Gd with large cross sections is the main response of the delayed 
signal in KASKA experiment, and we measure absorption probability and timing 
by using neutron sources. ^"^^Am-^Be source is useful for this purpose. In Am- 
Be sources, a particle produced by ^'^^Am decay is irradiated to ^Be nucleus, 
and then it forms excited state of ^■^C and decays to neutron plus 4.4 MeV 7. 
^^^Cf is a spontaneous fission source and can produce high rate neutrons. 

6.4 Calibration device manipulator 
Region-I 

A very basic calibration equipment is a vertical string suspendable in the LS. 
Calibration sources are put in the known position along the z-axis. This type 
of calibration will be most frequently performed. In KASKA experiment, three 
dimensional calibration is necessary to study the properties of signals near the 
acrylic walls. In order to accurately locate sources at target positions, we de- 
velop a calibration manipulator. We have two candidates of manipulator types, 
wire type like SNO and robot arm type. For the wire type, the area of the 
device that shadows the scintillation photons is small. However, for the wire 
type, the anchors hanging the wires may damage the acrylic vessel if the load to 
the anchor is too strong. In addition the anchors and wires must be kept in the 
LS. On the other hand, by the robot arm type, all components can be taken out 
of the LS for physics running, and the direction of sources can be controlled. In 
addition, by replacing the source with a CCD camera, we can perform visual 
inspections, during the oil filling and to measure absolute positions of the 3 
dimensional calibration device by comparing with the positions of LEDs which 
are distributed around acrylic sphere or PMT sphere. 



6.4.1 Wire type manipulator 

Fig.M7Ueft shows how to locate the calibration source in the detector by a wire 
type manipulator. The concept is to use a wire manipulator as in SNO detec- 
tor. Four ring-blocks are attached inside the inner acrylic sphere and a wire 
goes through the ring. One end of the wire is tied to the calibration device and 
the other end goes all the way up to the top of the chimney and the length are 
controlled by a stepping motor. By controlling the length of the three wires, 
the position in the horizontal plane is determined. The 4th wire directly comes 
down the chimney and is attached to the top of the calibration device. The 
weight of the calibration device is heavier than that of LS and by controlling 
the length of the 4th string, the vertical position of the calibration device is de- 
fined. The basic calibration is frequently performed to put a radioactive source 
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along the vertical line. Then less frequently 3 dimensional calibration will be 
performed. 



6.4.2 Robot arm type manipulator 

Fig.MTh'iglit shows how to locate the calibration source in the detector by a 
robot arm type manipulator. By using two joints and a rotating bar going 
through the chimney and a connecting arm, sources can be located at target 
points. In addition, by using the joint connecting a source to arm, the direc- 
tion of the source can be controlled. By making air chambers inside the arm, 
it is possible to make it zero-buoyancy to make the force to move the joint small. 




Figure 47: Source manipulator 



Region-II 

As for the Region-II, the main purpose of the calibration is to monitor the 
change of light output and attenuation length. String guides will be equipped 
in the Region-II as shown in Fig.-BSl A string runs through the string guides. 
Radioactive sources are connected to the string and moved along the vertical 
direction by controlling the length of the strings. 

Region-Ill 

Low luminous LEDs, high luminous infrared LEDs and CCD cameras are mounted 
on the inner surface of the stainless steel tank. The low luminous LEDs are used 
to provide one photon level signals to PMTs to monitor the PMT gain easily. 
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High luminous LEDs are placed at known positions and become a fiducial marks 
for the CCD cameras. The CCD cameras are used to watch the condition of 
the detector and to measure the position of the calibration sources making use 
of LEDs as reference positions. 




Figure 48: Region II calibration 



7 Backgrounds 

The background reduction and its precise estimation is a main issue of KASKA 
experiment because the uncertainty associated with background estimation is 
one of the most severe sources of the systematic error and the detector system 
has to be designed carefully to eliminate it. 

FigEHl shows an example of the background measured in CHOOZ experi- 
ment. 

The events in region-D are caused by accidental coincidence of 7-rays from 
detector material and rocks. Those in Region-C are interpreted as correlated 
background caused by fast neutron. The events in region-A are caused by 
accidental coincidence between 7-rays and high energy signals such as recoil 
protons. The events in region-B can be associated with decaying muon after 
stopping in the detector. Those backgrounds leak in the neutrino candidate 
region and CHOOZ experiment observed 4.0% of correlated background and 
1.7% of accidental backgrounds. 

The Kashiwazaki-Kariwa nuclear power station has 7 reactors and it is im- 
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Reactor CTJ 



Reactor OFF 




Figure 49: Correlation of prompt energy and delayed energy measured by 
CHOOZ group Pni- The left data were taken during reactor ON and the right 
data, during reactor OFF, showing backgrounds. The events in Region-A and D 
are caused by accidental coincidences, Region-B are caused by a Michel electrons 
from stopping /x and Region-C are supposed to be caused by the fast neutron 
background. 



possible to take reactor-off data. The strategy of the KASKA experiment for the 
background evaluation is to reduce the absolute background rate to be less than 
a few% level and then to make the absolute error of the background estimation 
smaller than a half %, by making use of the following methods. 

(1) Powerful reactors: The Kashiwazaki-Kariwa nuclear power station is the 
world most powerful reactor complex whose thermal power is 24.3GW, which 
produces 3 times more neutrinos than CHOOZ power station. With this high 
neutrino flux, the S/N ratio and statistics become large at a given detector 
configuration. 

(2) Underground laboratory: The neutrino detectors are placed in a tunnel 
excavated at the bottom of deep shaft holes to reduce the cosmic-ray rate. 
The most severe backgrounds, such as fast neutron and correlated signals from 
spallation products are originated by cosmic rays. The water- equivalent depth 
of the far detector is 260 mwe, which reduces the cosmic ray rate down to 1/240 
of that on the surface. For near detectors, the depth is roughly 90 mwe which 
reduces the cosmic-ray rate down to ~l/30. 

(3) Gd loaded liquid scintillator: The total energy of the 7-rays from the 
neutron captured Gd is ^ 8 MeV. This is much larger than the highest natural 
7-ray energies of 2.6 MeV and natural highest /3-decay energy of 5 MeV. Thus 
the delayed signal is free from background of natural radio isotopes. 

(4) Thick shields: The external backgrounds such as 7-rays and fast neutrons 
can be reduced by thick shields. Fig.-ISOl compares the KASKA detector and 
the CHOOZ detector. The thickness of the KASKA shield is much thicker than 
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that of CHOOZ detector and the external backgrounds are significantly reduced 
compared with the CHOOZ case. An option is to make the buffer oil as a weak 
scintillator having a light output of 1/10 of Gd-LS. 90 cm thick buffer between 




Figure 50: Comparison of the CHOOZ and KASKA detector. 

PMT array and the liquid scintillator reduces the 7-rays from PMT glass. There 
is at least 2.2m thick buffer and water barriers between the liquid scintillator 
and iron tank. The 7-ray and neutron background from outside of the detector 
is reduced to negligible level with these buffers. 

(5) Cosmic-ray tracker: When cosmic-ray hit ^^C, radioactive elements are 
produced by spallation reactions. Among them, long lived /?-neutron emitters, 
such as ^Li or ^He are most significant background sources. Because spalla- 
tion products are produced along the cosmic-ray track, the rate of such kind 
of background can be estimated by the distribution of distances between the 
cosmic-ray track and signal candidate. A cosmic-ray tracker system is imple- 
mented to estimate this kind of spallation backgrounds. 

7.1 Background Sources 

Background signals are caused by radioactive elements such as ^^^U(t = 4.5 x 
10\), 232Th(r = 1.4 X 10^°y), ^°K{t = 1.3 x lO^y) and 6"Co(r = 7.6?;). ^asy, 
and '^"K naturally exist in various materials due to their extremely long 
lifetime while ^°Co is artificially contaminated in steels in smelting process. 
Fig. 1511 shows the decay chain of the ^ssy g^^^ 232rpj^^ j£ radioactive elements 
decay within liquid scintillator, both f3 and 7 rays contribute as background 
(Internal background). The a particles can not be a background because the 



60 



i «'l.l82MeV 

■!94Th(2M)->3:wPaa.2nJ^=:*-!L'l?.r,EE) 

p eaSl $2.20 i a 4.768 

J, a 4.6Sr 
S6Ra<t8£i0yJ 
la -1,777 

I u 5 -182 

J,a600i 

■**R>g!»0-»=''BLi,?0[ji)^-"Paio>J 
pi.03 p3.2G iaZeSO 

.■i."Pbe2yi^-"'B](i:-.i!..ij^='"Pul!,«d) 

PO-P0i I o, 52.98 

«*Tha.4El03r) 

I a 4.007MeV 

«2«Sa(s.Bj-)->-'Ac(ti,iw-»=="Th{i.!)y) 
pp.046 P3.127 la 5.420 



i a 6.282 

la6.V7 |B,2iC6SS 

pO.S73 J a6 05(34%) la8.';SS 
a*n6.1in)-»2'«Pb.. 



Figure 51: Decay of ^^^U chain, ^^^Th chain and *°K. 
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visible energy of the a's produced in natural radio isotopes are quenched to less 
than 0.5 MeV. If these elements decays outside the liquid scintillator, only 7- 
rays contribute as background since the range of a few MeV /3 is only a few cm 
(External background). Table ITTI shows their background rates which satisfy 
the positron energy window. The highest decay energy of the /3 decay is 5 MeV 
(^°^T1 ^208 p|-, 232rp-|^ chain) and these radioactive elements do not satisfy 
the neutron signal criterion of E > 5 MeV. The backgrounds for Gd signal are 
caused mainly by actual neutron absorption, where the neutrons are produced 
by cosmic-ray associated reaction with materials. 



Table 14: U, Th, K, Co activities 



Element 


Bq/g 


#7/decay 
(E>0.7MeV) 


#7/g 


# (7 + /3)/docay 
(E>0.7MgV) 


# (7 + /3)/g 


238 u 


12,400 


2 


25,000 


3 


37,000 


232Th 


4,100 


1.5 


6,200 


2 


8,200 


40k 


255,000 


0.11 


28,000 


0.5 


130,000 






2 




1 





Table 15: Internal background 



Component 


Quantity 
(kg) 


238 u 
(PPt) 


232 rp]^ 

(ppt) 


40k 

(ppt) 


/3 + 7 rate 
(Hz) 


Scintillator Oils 


25,000 


< 0.1 


< 0.1 


< 0.1 


< 0.5 


PPO 


150 


4 


< 50 


4.2 


< 0.17 


bis-MSB 


2.5 


48 


60 


6.8 


0.008 


Gd 


6.1 


8.4 


12 


< 1200 


< 0.96 


Acrylic 


700 


8 


50 


8.2 


1.2 


Internal Total 










< 2.9 



7.1.1 Internal Background 

* The liquid scintillator 

The amount of oil used for both LS-I and LS-II is 25 tons. The liquid scintil- 
lator will be formulated by dissolving PPO, bis-MSB and organic Gd into an 
oil which is made mainly of pseudocumcne and paraffin oils. The U/Th/K con- 
tamination of these materials has been measured by ICP-MS and flame atomic 
absorption method. The oils are produced from a number of steps of distillation 
and refinery processes and are intrinsically pure. The U/Th/K contaminations 
are less than 10~^^g/g and the total single rate from oil is less than 0.5 Hz. The 
single rate from PPO, bis-MSB and Gd samples were estimated to be less than 
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1.2 Hz from its measured contaminants. 
* Acrylic vessels 

The weights of the inner and outer acrylic vessels are 240 kg and 900 kg, re- 
spectively. Because half of the 7-rays from the outer acrylic vessel do not come 
in the scintillator region, the effective mass of the outer vessel is 450 kg and the 
effective total acrylic weight is about 700 kg. The single rate from acrylic tanks 
are estimated to be 1.2 Hz. 

In summary, the internal background is expected to be < 2.9 Hz as shown in 
table O 

7.1.2 External Background 

For external background, only 7-rays can enter liquid scintillator. 



Table 16: external background |W 



Material 


Quantity 
(ton) 


238 u 
(ppb) 


232Th 
(ppb) 


40k 
(ppb) 


60 Co 
(mBq/kg) 


Hit rate 
(>0.7MeV)(Hz) 


PMT Glass 


0.3 


46 


120 


10 




1.2 


PMT support 


2 


1.9 


4.7 


0.08 


9 


0.2 


S. S. Tank 


32 


1.9 


4.7 


0.08 


9 


0.05 


Iron Tank 


300 


1.9 


4.7 


0.08 


9 


0.01 


Soil 




3,000 


8,000 


2,000 




0.15 


Concrete 


180 


2,800 


5,400 


1,700 




0.43 


Total 












2.0 



* PMT Glass 

In KASKA, 300 10" PMT will be used. The total weight of the PMT glass 
is 300 kg. For Hamamatsu low background glass, U/Th/'"'K contaminations 
are measured to be 46 ppb, 120 ppb and 10 ppb, respectively and 650/s of 7- 
rays which satisfy the positron energy cut are generated from the glass. After 
penetrating 70 cm thick buffer oil, the single rate in the LS becomes 1.2 Hz. 

* 7-rays from the stainless steel and iron structure. 

The mass of the stainless steel tank, PMT and steel tank is 30 tons, 2 tons and 
300 tons, respectively. Generally a steel contains lOmBq/kg of ^"^Co as well as 
U/Th/K. ^"Co emits two 7-rays whose energies are, 1.17 MeV and 1.33 MeV. 
After passing through at least 90 cm of buffer oil, the single rate in the LS 
becomes 0.25 Hz. 

* gamma rays from soil & concrete: 

The radioactive contaminations in the soil at the near-B point are measured in 
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a boring study to be, U/Th/K=3/8/2 ppm. A simulation shows the single rate 
of 7-ray hit from soils and concrete is only 0.6 Hz 

The total single rate which is originated from natural radio activities is ex- 
pected to be < 5 Hz. Fig.l52lshows the simulated external background spectrum. 




Figure 52: External background spectrum ^S] 



7.1.3 Cosmic Rays 

In order to reduce cosmic-ray associated background, the detectors are placed 
at 150m depth underground for far detectors and at 50m depth for near detec- 
tors. The mean density of the soil around Kashiwazaki-Kariwa Nuclear Power 
station is 1.75g/cm'^. 30 weight % is water component. The water equivalent 
overburden at 150 m underground is 260 meter- water-equivalent (mwe) and that 
at 50 m undergaround is 90 mwe. Fig-ISll shows the cosmic ray flux at various 
depth calculated by an empirical formula j48| : 

I{h) = 1.4 X ^ih + 11)^'-^^ exp(-7.0 x 10-^/i)[/mVs] (18) 

where h is the depth measured in mwe and the factor 1.4 is a normaliza- 
tion factor which is chosen to match to our simulation result at 150m depth. 
Fig. 1541 shows the energy spectrum and zenith angle distribution at far detector 
depth At 150m depth, the cosmic ray flux is ~ 0.55/m^/s and at 50m depth 
it is ~ 6.6/m^/s. The cosmic-ray hit rate at inner detector region (region-I, II, 
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Figure 53: Cosmic Ray Rate 




Figure 54: Cosmi ray energy spectrum (left) and zenith angle distribution (right) 
at far detector depth [17|. 
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Ill) is ~15 Hz for far detectors and ~170 Hz for near detectors. The hit rate at 
outer detector region (region-IV, V) is ~40 Hz for far detectors and ~500 Hz 
for near detectors. The mean energies of the cosmic rays are ~70 GeV and 
^^20 GeV, respectively. 

The cosmic-ray veto system consists of at least 1.6 m thick buffer oil in 
region-IV for vertical direction and 0.7 m thick buffer oil in region-IV and 0.5 m 
thick water in region-V in horizontal direction. A minimum ionizing particle 
deposits at least 260 MeV of energy in region-IV at vertical direction, 110 MeV 
of energy in region-IV and 100 MeV of energy in region-V at horizontal direc- 
tion. Together with the water cherenkov veto and cosmic-ray tracker, >99.5% 
of cosmic-ray rejection efficiency is expected. A fast neutron which collides with 
a proton in the region-IV and deposit visible energy of few MeV is also expected 
to be vetoed by cosmic-ray anti counter. In order to avoid the instantaneous 
cosmic-ray associated backgrounds, 0.5 ms of dead time will be applied after 
each muon passage in the liquid scintillator region and 0.2 ms of dead time is 
set after the muon passage in non scintillating region. The livetime inefficiency 
introduced by those dead time is around 1.6 % for far detectors and 20 % for 
near detectors. The Ve/ 1^- rate is higher for near detectors than far detectors and 
the mean muon energy at near detector depth is only 1/3 for far detector depth. 
Therefore the N/S of cosmic-ray associated background in near detectors is a 
half of that of far detectors. 

The liquid scintillators consist of hydrocarbonic oils and chemicals which 
consist of hydrogen, carbon, nitrogen and oxygen. When a cosmic ray hits 
such nucleus, it may break the nuclei producing radioactive nucleus. However, 
the atomic components of the LS's are roughly H:C:N:0=1.8 : 1 : 0(10^'^) : 
O(10~^) and carbon contributes mainly to such spallations. An example of 
such reaction is, 

/i -H^^ C -> -f " C* -f n (19) 

The excitation of the produced nucleus decays instantly and its signal is killed 
by cosmic-ray veto. However, ^^C decays with electron capture with half life 
of 20.4 ms, which is long enough to survive the cosmic-ray veto. The neutron 
produced will be absorbed by proton or Gd, which produces a delayed signal. 
Fig-Eilshows a list of long lived radioactive elements which can be produced by 
spallation of carbon nucleus. Table [T7I shows their production rate in region I 
and II, calculated by using the cross sections taken from and assuming that 
the cross section is proportional to E^{Ge\^Y'^'^^ , where is a mean energy of 
cosmic-rays. The total production rate within region-I and II is 0.036/s for far 
detector and 0.27/s for near detectors, both are much smaller than 7-ray rate 
from detector materials. If the neutrons produced by the spallation reaction is 
absorbed by Gd, they generate delayed signal. The rate of such neutron signal 
which survives the cosmic-ray veto time is, 



66 








H 

1 


He 


Li 


Be 
4 


B 
5 


C 

6 







H(99.99) 








stable 


5nii5T?ibl(? 


1 


n 

iO.eim 

b-(0.7S2) 


2H(0.[il) 


3He(l.JE-6) 






(abundance) 


name 
Lifetime 
mode 


2 




511 

12.3y 


—mm 










3 










7S^ 

53.3d 
EC(0.S62) 


EB 

0.77s 
EC2a{1.798) 


?c 

0.127i 
ECp2a(16.50) 


4 






6He 
0.808s 


7Li(92.5S) 


SBe 
0.0G7fs 

:a 


9B 
0.855 


IOC 
15.1SS 
EC(.3&18) 


5 








0.844s 
b-2a(1.60) 


?Be{iflCil 




lie 

lO.iSras 
EC{1.98) 


6 






SHe 
0.122s 
6-11(10.65) 


0.17S5 
b-ii(lj.61) 


lf)Be 

l,6E6y 
b-{0.5511 






7 










llEe 
13.Ss 
b-a( 11.51) 


o.o:u3 

b-ja{13.37) 





Figure 55: Long life radio isotopes possibly produced by cosmic-ray spallation 
reaction of carbon |49|. 



Table 17: Production rate (Region-I+II) 





Far 


Near 


mean energy (GeV) 


68 


19 


/I flux (/m^/s) 


0.55 


6.6 


"C{20.4ms) 


2.6 X 10-2 


1.9 X 10-1 


^Be(53.3d) 


5.8 X 10-^ 


4.8 X 10-2 


"Be{13.8s ) 


< 5.5 X 10-5 


< 5.0 X 10-* 


iOC(19.15s) 


3.5 X 10-^ 


2.4 X 10-2 


8Li(0.844s ) 


1.4 X 10-* 


8.4 X 10-* 


''He{0.808s) 


4.6 X 10-* 


3.3 X 10-3 


*B{ 0.77s) 


1.9 X 10-* 


1.5 X 10-3 


9C(0.127s) 


1.6 X 10-* 


7.5 X 10-* 


^Li(0.178s)+*He(0.122s) 


6.0 X 10-5 


2.8 X 10-* 


sum 


3.6 X 10-2 /s 


2.7 X IQ-i /s 
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where Nn is the number of neutrons produced in the spallation reaction, // is 
production rate of radioactive spallation nucleus in region-I and T^^io is the 
cosmic-ray veto time. (3 is the fraction that a neutron produced in region-II 
'spills-in' to the region-I and is absorbed by Gd, which is an order of 10%. For 
the case of ^^C and T^q^q = 0.5 ms, /„ < 3.2 /day for far detector. For other 
nuclei, it is much less than that and total neutron rate is less than 5 /day. 

As described in the next section, there is a case that fast neutron produced 
by hard collision of /.t and nucleus will be thermalized in the liquid scintillator 
and eventually absorbed by Gd. According to a MonteCarlo simulation, this 
kind of background takes place at a rate of 2.2/day in far detector assuming 
99.5% efficiency for muon veto. 

7.2 Accidental Backgrounds 

Accidental background is the events in which independent background signals 
accidentally happens within the delayed coincidence timing window and re- 
garded as neutrino signals. 

The accidental background rate Jab is expressed as, 

lAB=fe+fn^ (21) 

where and /„ are rates of e+ like signals and neutron like signals, to is the 
coincidence window. For KASKA experiment, the energy windows for positron 
and neutron signals are, 0.7 MeV to 9 MeV and 5 MeV to 11 MeV, respectively. 
The timing window is 200 /iS. 

Table 1181 summarize the single rate and accidental background rate for far 
detector. Accidental background rate can be estimated precisely by opening 
coincidence window randomly. 

Table 18: Single rate and accidental background summary. 



Source 




fn (/day) 


Internal 


<2.9 




External 


2.0 




Cosmic ray spallation 


0.036 


< 7.2 


sum 


< 5 


< 7.2 


Accidental background 


<0.01/day 



For near detectors, S/N is expected to be better than far detector because 
the mean energy of cosmic-ray is much smaller and the reaction cross section is 
smaller. 

7.3 Correlated Backgrounds 

Correlated background is a signal in which a single event produces two contigu- 
ous signals, each satisfies the prompt and delayed energy windows for De selec- 
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tion within coincidence time window. Fast neutron background and /3-neutron 
emitting radio isotopes are categorized as correlated backgrounds. 

7.3.1 Fast Neutron Background 

When a fast neutron collides with proton in a scintillator, there is a chance 
that the visible energy of the recoiled proton falls in the energy window for 
prompt signals. Also there is a chance that the original neutron looses its energy 
scattered by protons around and eventually absorbed by Gd. In this case the 
fast neutron signals are misjudged as Ue events. Such fast neutrons are often 
produced by inelastic scattering of cosmic-ray muon, 

H + H + n + X (22) 

According to the KASKA simulation, the fast neutron background (FNB) 
coming from the muons which go through outside of the veto system is negligi- 
bly small thanks to the thick buffer regions. Most of the FNB comes from the 
muons which come in the detector region but can not be detected due to the 
inefficiency of the veto system. For 99.5% of veto efficiency, 0.2 /day after the 
event selection, or ~ 0.7% of the Vf, signal is expected from the simulation. The 
rate of FNB can be estimated by number of ways. Because most of the FNB is 
associated with the veto inefficiency, the rate can be estimated by applying veto 
inefficiency artificially when analyzing the data. It can be estimated from the 
energy distribution of the prompt signal in the non-reactor neutrino energies or 
spallation background energies, such as beyond 15 MeV and less than 1 MeV. 
Fig. [S^ shows visible energy distributions calculated by a MonteCarlo simula- 
tion. The figure shows the energy distribution is flat down to 4 MeV and there 





prompt signal 


■ .... 1 1 .... 1 


Entries 2847 
Mean 10.57 
RIVIS 8.632 







5 10 15 20 25 30 

prompt deposit (MeV^ 



Figure 56: Visible energy distribution of delayed energy |47| . 

is an enhancement at i? < 4 MeV. The enhancement can be estimated from 
the event rate below 1 MeV. It will be possible to understand the energy distri- 
bution with real data by tagging cosmic-rays. The radial position dependence 
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of the events also gives information about fast neutron background since the 
neutron backgrounds tend to occur near the detector surface. The Palo- Verde 
type liquid scintillator and 7 catcher liquid scintillator have nentron/7 pulse 
shape discrimination capability and this feature may be of help to separate the 
neutron and positron signals. Assuming the rate of fast neutron background 
can be estimated with <30 % accuracy, the absolute error of this background 
estimation is <0.2 %. 

7.3.2 Spallation Background 

Some spallation products have decay mode to produce P + n. In such decays, 
the P signal mimics the prompt signal of the reactor neutrino and the neutron is 
absorbed by Gd producing exactly same signal as reactor neutrino delayed one. 
Among these P + n emitter, ^Li and ^He are most serious background source 
because they can be produced from ^^C and have long lifetimes and can not be 
rejected by the veto time applied after muon passage. 

KamLAND group showed that ^Li dominates and contribution from *He 
was less than 15% in their data taken at 1000m deep underground and its 
amount is consistent with reference- '50' . According to the reference, the pro- 
duction rate of ^Li-|-*He in region-I is 1.2 /day for far detectors. Because the 
branching ratio of the f3 — n decay mode is 50% for ^Li and 16% for ^He, the 
neutron associated background rate becomes at least half, which corresponds to 
<1.5 % of the reactor i/g events. The other (3 — n emitters are assumed to be 
at least O(10~^) times smaller because they have to be produced from oxygen 
nuclei. The ^Li background rate can be estimated with following methods. 
*Energy Spectrum: Because the Q value for neutron associated decay of ^Li is 
11.2 MeV, this background can be estimated from the number of events between 
9 MeV to 11 MeV. ~5% of ®Li signals falls in this non energy window. 
* Cosmic-ray track: Because decay mode of the ^Li is (3 decay, the prompt like 
signal occurs within a few cm along the cosmic-ray track. 

The production cross section and its extrapolation on mean muon energy was 
measured with accuracy 20% [SU]. Assuming the background rate can be esti- 
mated with precision <30%, the absolute error comes from this background is 
<0.5%. 



7.3.3 Stopping /i 

Fig. 123 shows the probability for cosmic-rays to stop after passing through 
1 g/cm^ of material. In order for cosmic rays to satisfy visible energy criteria 
of prompt energy selection, the muons have to stop at less than 40 cm from the 
PMT layer if weak scintillator is used for the buffer oil. The probability to stop 
less than 40 cm in the region-Ill oil is 0.4% and absolute rate is ^200/hr for 
far detectors. Assuming 99.5 % of cosmic-ray veto efficiency, the survival rate 
becomes ^1/hr. However, the stopping position is still in region-Ill and can 
easily be distinguished from i/g events which occur within region-I. 
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Figure 57: Cosmic-ray stopping probability. 



Most of stopping muons decay as, 

/X — > e + + j/e (23) 

The maximum kinetic energy of electron or positron is 53 MeV which produce 
^^5 MeV of visible energy in region-Ill. Some of these events satisfy the delayed 
signal criteria. Assuming 10% of delayed energy cut efficiency for the Michel 
electron, and timing cut efficiency of e^^^^/^-^'^*, the correlated event rate is 
1.5 event/day. For fj,~ case, 7.9% of stopped fi" is captured on -^^C and 80% of 
which goes to unstable nucleus which typically emits neutron |14j . 

fi- +^^C^n + X + z/^ (24) 

If the stopping muon, which satisfies the prompt event selection, emits neu- 
tron, it can be a correlated background. Assuming 1 % of neutron produced 
by such stopping muon is captured by Gd, the rate of this background is, ~2.5 
events / day. 

However, these events can be easily distinguished from z/g events making use 
of the extreme position difference of the signals and reduced to a negligible level. 

Table IT^ summarizes the expected background rate and its estimation error. 
For near detector case, the S/N is two times better than the far detector 
case and the absolute error after near/far comparison is 0.6% or better. 



8 Systematic Error 

It is difficult to accurately estimate systematic error before experiment starts. 
In this section an upper limit of KASKA systematic error is estimated by com- 
paring the recorded errors by CHOOZ experiment. 
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Table 19: background summary 



source 


absolute 


error 


Accidental 
Fast Neutron 
9Li+«He 
Stopping /X 


<0.02% 
0.7% 
<1.5% 
nagligible 


negligible 
<0.2% 
<0.5% 

negligible 


Combined 


<2.2% 


<0.5% 



8.1 Near/Far cancellation 

If neutrino oscillation takes place, the expected neutrino spectrum; fJ-uiE) at 
distance L from a reactor is calculated as 

fi,{E, L) = fi°iE, L){1 - sin^ 2^13 sin^ ^^^iE, L)) (25) 

where ^1(E,L) is the energy spectrum in case of no neutrino oscillation calcu- 
lated by the equation below. 

t^"(E,L) = ^Npa.p{E)ME) (26) 

where Np is the number of target protons, fiy{E) is the number of generated 
neutrinos in reactors and a^p{E) is the differential cross-section of the inverse 
neutron /3 decay, 

+ P e+ + n (27) 

In the KASAK experiment, the oscillation will be searched for mainly by using 
the deficit of the total event rate. For such case, the expected number of events; 
Ail. is, 

/i, =/z^(l-sin2 20i3A(L)) (28) 

where, 

^,l = J iil{E,L)dE (29) 

is the total number of events in case of no neutrino oscillation. 

t^l{E,L)s\n'<i>,:,{E,L)dE 

is the deficit of event at distance L from the reactor. Fig.l58lshows A{L) for the 
case of Am^ = 2.4 x 10^'^eV^. At the first oscillation maximum of L ^ 1.8km, 
only 20% of De remains. 

For simplicity, we assume a case to measure sin^ 2^13 with a scheme of one 
reactor, and one far detector at a distance L, one near detector at a distance 
L'. If we observe Nobs events in the far detector and N^f^^ events in the near 
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Figure 58: Deficit probability (A(L)) 

detector, we consider e/.t^ = Noh^ and e'^i'^ = N'^^^, where e and e' represent 
detection efficiencies of far and near detectors, respectively. Therefore, the most 
probable sin^ 20i3 is. 



where 



(32) 



is expected to unity in case of no neutrino oscillation. It is noted that a^p and 
fi, do not appear in Ea. H31() and the measurement accuracy does not directly 
depend on the errors on such parameters. If value of Ea. H31|l is significantly 
larger than 0, it indicates that non-0 sin^ 2^13 is observed. 

Since sin^ 2^13 is known to be small, the measurement error is, 

where AA = A(i) - rA{L') sa 0.7 for KASKA base hues. In order to obtain 
a small error, having a large AA is essential. Sunmrr means uncorrelated error 
for both detectors. By taking the ratio, the common uncertainty components 
cancel out and only uncorrelated components remain. The term of square root 
is the statistic error. 

In the KASKA experiment, the main event selections are, 

(1) 0.7 MeV < Eprompt < 9MeV 
(2) 5MeV < Edeiay < llMeV 

(3) ^^ls <At< 200ns 

The detection efficiency is multiplication of each event selection. 

e = e^+CneAt (34) 
Errors come in from the relative uncertainties of each parameter. 
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8.2 Detector Associated Error; 5, 



uncorr 



{eN,)/eN, 



8.2.1 Number of target protons; Nj 



p 



The number of target proton is calculated by, 



Np = ModLSPp 



(35) 



where pp is the proton density per unit LS mass and McdLS is the mass of the 
Gd-LS within the fiducial region. The Gd-LS for all the detectors will be stored 
in a storage tank before introducing in the detectors and sent to each detector 
from the tank. So that in principle, the component of the LS is same for all the 
detectors. 

KASKA does not apply fiducial cut because it introduces a large uncertainty. 
Instead, the reactor neutrino signal is identified by existence of the neutron- 
captured Gd signal, regardless where the prompt signal occurs. Because region- 
II scintillator does not contain Gd, there is a natural cut-off of the fiducial region. 
The main error of the neutrino detection efficiency comes from the relative 
difference of the fiducial mass of the Gd-LS between the near/far detectors. 
The mass will not be defined by the volume of the container because it possibly 
deforms under the pressure of the liquids. Instead, it will be defined by total 
mass of the LS which is measured precisely when filled in the detector. The 
Gd-LS will be delivered by 200 liter oil drums and the weight difference before 
and after the filling is measured with 0.3% accuracy weighing machine. The 
three weighing machines used in the far and near detectors are cross calibrated 
occasionally during the filling process. The diameter of the chimney is 30 cm 
and 0.1% of volume error corresponds to 10 cm of head level difference, which 
can easily be detectable by a level monitor. 

If the neutron produced in a D,, event goes out from region-I, it introduces an 
inefficiency (spill-out effect). Contrarily even if a event occurs in region-II, 
it is considered to be Pe signal if the produced neutron comes in the region-I 
and is absorbed by Gd (spill- in effect). The spill- in effect largely cancels the 
spill-out inefficiency. Fig. -1551 depicts such effects. The error introduced due to 
the edge effects is 



where ai,aii and pi,pii are diffusion length of neutrons and proton density 
in region-I and II, respectively. R is average radius of region-I. The main error 
component comes from uncorrelated uncertainty of the diffusion length. Be- 
cause diffusion length is insensitive to liquid scintillator components, the error 
is assumed to be negligibly small. 

In summary the relative uncertainty of the fiducial mass is safely claimed to 
be better than 0.5%. 



AN, 



p 



3(Acrj - {pii/pi)Aajj) 
2R 



(36) 
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Figure 59: Radial dependence of detection efficiency. Average radius of 1.1 m 
is used as a typical size. Diffusion length of neutron is assumed to be 6 cm. 



8.2.2 e+ Energy Cut Efficiency; 

The threshold for the positron energy is 0.7 MeV which is Scr below the positron 
annihilation energy of 1.022 MeV. Fig. 1601 shows the reactor neutrino spectrum 
near the energy threshold. Even if the energy resolution is cte/E = 15% the 
rate below 0.7 MeV is less than 0.03% and even if 0.1 MeV of detector-relative 
systematic bias takes place, the relative difference in the number of event is less 
than 0.1%. 
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Figure 60: Energy spectra near the threshold. The dashed line shows original 
energy spectrum and the solid like shows the spectrum smeared by finite energy 
resolution of ue/E = 15%. 
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8.2.3 Neutron Capture Efficiency; e„ 

The efficiency for neutron signal cut is a multiply of neutron absorption efficiency 
on Gd; tabs and the detection efficiency of Gd 7 signal; ecd-y , 

En — ^abs^Gdj- (37) 

The neutron capture efficiency for Gd is expressed as, 

eod = (38) 

PGd cr^a 

where p is the number density of the corresponding element, cr" is the 
neutron absorption cross section as shown in the table 1201 
The ratio of the absorption cross section is (^abs/'^abs — ^-^^ x 10^^. 



Table 20: neutron absorption cross section 



Element 


Q value (MeV) 


Abundance(%) 


cross section (b) 




2.2 


99.985 


0.3326 ± 0.0007 




8 


14.80 ±0.05 


60900 ± 500 




8 


15.65 ±0.03 


254000 ± 800 


Gd Average 






48800 ± 400 



The atomic component of the KASKA liquid scintillator is 

#H:#C:#0:#Gd=l. 85:1:0.018:0.00009 

and p^ts/Pabs — 2-01 x 10*. Therefore, the neutron absorption cross section 
on Gd is, eod ~ 0.88. The capture efficiency error associated with the Gd 
concentration error is, 

^ = (1 _ ecd)^ = 0.12^ (39) 

EGd Pad PGd 

Therefore the capture efficiency is insensitive to the change of the component 
of the liquid scintillator. Because the Gd-LS for all detectors is taken from the 
same storage tank, the variation of the component is supposed to be small and 
the error on ecd is estimated to be better than 0.5%. 



8.2.4 n Energy Cut Efficiency; e„ 

The total energy of the gamma rays from neutron captured Gd is 8 MeV. The 
energy threshold for the neutron signal is 5 MeV, which is low enough and the 
uncertainty of the threshold does not cause error of the event selection. Fig. ICT 
shows the Gd-captured neutron signal observed by the CHOOZ group (10) . 

CHOOZ group assigned 0.4% as absolute error with energy threshold 6 MeV. 
Because the energy threshold can be set lower than CHOOZ case, thanks to a 
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Figure 61: Gd-captured neutron signal observed by the CHOOZ experiment [THj 

new buffer layer in front of PMT array, and energy containment is better, and 
because error of KASKA is relative error, it is expected to be much better than 
0.4%. The threshold energy can be calibrated by n +^'^ C -^^^ C + ^{bMeV) 
events. 



8.2.5 Efficiency for the timing cut; ea* 

The efficiency for the timing cut is. 



1 _t2 

eAt = - e-di = e- -e- (40) 

T .It, 



and its error is, 



Se 1 ,tie^~ — t2e^~ At _tx Ati _^i2At2. 
— = — Tl —i 1 + e - — ) (41) 



e T — e 



T T T T 



For KASKA, r w 30/is, ti = 1/^s and t2 = 200^sec. So that. 



EA, = 0.966 (42) 

and 

^^ 0.058^-1.0^ + 0.0014^ (43) 
err r 

The absolute timing will be quantized in unit of system clock frequency 
such as O(lOns) and the timing of scintillation signal is O(lOns). Therefore the 
error induced by absolute timing is negligibly small, 0(0.03%). Because the 
components of the Gd-LS is same for all the detectors, essentially At is small. 
The absorption time r can be measured in — situ with real data with precision 
of a few % and the At term is order of 0.1%. 
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8.3 Error from Reactor Neutrino Flux; f^, 

The reactor neutrino flux; f^, is known to accuracy of 3.4%. For one reactor and 
2 detectors case, tfie uncertainty completely cancels. However, for KASKA case, 
there are 7 reactors and three detector locations whose near/far baseline ratios 
are different. In this case the cancellation becomes imperfect and some ambi- 
guity remains. This effect is estimated analytically in ref. The following is 
an intuitive description. The neutrino flux is calculated as, 

fu ^ Pth'x nf X r^^ (44) 

where Pth is the reactor power, Uf is the fission rate per thermal power and 
rjy is the neutrino flux per flssion. Their uncertainties are, 2%, 1%, and 2.5% 
respectively. However, 77^ is common to all reactors and cancel out by near/far 
comparison. All the reactors in Kashiwazaki-Kariwa nuclear power station will 
be more than 10 years old and their fuel components have became equilibrium . 
Thus the ratio of the fission rate of each fissile element are similar when averaged 
over a few years of run period. Although the original error is ^^1%, it will be 
canceled out to the level of negligibly small. The ambiguity from the baseline 
difference comes from the error of thermal power measurement, multiplied by 
the standard deviation of the square of the ratio of far/near baselines. 



M P\^^\ R 



where Rr is the square of the ratio of far and near basehnes, Rr ~ {LrjL'y)'^ 
and R is the average of Rr\ 

R^\Y.^- (46) 

r 

For KASKA case, the value of the square root in Ea. l|^ ~ 0.07 and the 
error introduced from the neutrino flux uncertainty is estimated to be 0.14%. 



8.4 systematics summary 

Since the current upper limit of sin^ 2^13 was measured by CHOOZ experiment 
and the types of systematic uncertainties of KASKA experiment are similar 
to those of CHOOZ, it is appropriate to compare systematic uncertainties of 
both experiment to see the improvement of the sensitivity. Table [2] and 1^ 
compares the systematic errores so far discussed. The total systematic error 
can be controlled to be better than 1% and possibly be as low as 0.5%. The 
sin^ 26*13 sensitivity shown in Fig.-|221are calculated with these estimations. 



9 Conclusions 

The KASKA project enjoys world's most intense low energy v^. flux from Kashiwazaki- 
Kariwa nuclear power station and has opportunity to pin down several neutrino 
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Tabic 21: Detcx-tor Associated Systoinatics Errors. 





CilOOZ 


KASKA 


e+ energy {e^+ ) 


0.8% 


<0.1% 


e"*" position 


0.1% 




n capture (eabs) 


1.0% 


<0.5% 


n energy (ecd-,) 


0.4% 


<0.4% 


n position 


0.4% 




n delay (cAt) 


0.4% 


<0.1% 


e+ — n distance 


0.3% 




n multiplicity 


0.5% 






0.8% 


<0.5% 


Combined 


1.76% 


<0.8% 



Table 22: Systematics Summary 



source 


CHOOZ 


KASKA 




1.9% 


negligible 


Detection efficiency (e) 


1.76% 


<0.8% 


Reactor Power (Pth) 


0.7% 


0.2% 


Energy released per fission 


0.6% 


negligible 


Background 


0% 


<0.6% 


Combined 


2.7% 


<1.0% 
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oscillation parameters, such as ^13, precise 612 and Am^^. Moreover, by com- 
bining with accelerator-based neutrino oscillation measurements, there are pos- 
sibilities to find non-0 sin 5 and to solve 623 degeneracy. 

The first phase of the project is to measure sin^ ^13 with sensitivity ten times 
better than the current limit. The detector R&D has progressed well and it will 
be possible to start data taking 3 years after start of construction. 
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